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I.  INTRODUCTION  AND  SUMMARY 

The  prime  event  during  this  Interval  was  the  conference  In  London  as 
sponsored  by  the  Rank  Foundation.  As  mentioned  In  the  last  quarterly.  It  was 
devoted  entirely  to  the  subject  of  scanning  microscopy.  Our  work  was  given  a 
great  deal  of  credit  by  the  various  speakers  --  all  of  whom  recognize  some  of 
the  Inherent  advantages  of  scanning  for  the  various  forms  of  microscopy  — 
acoustic,  optical,  and  x-ray.  The  electron  microscope  was  not  Included  since 
It  hs  widespread  coverage  elsewhere. 

The  book  containing  the  papers  has  now  appeared.^  It  has  two  Important 
features:  1)  It  contains  a  good  deal  of  valuable  reference  material  on  this 
field,  and  2)  It  was  printed  within  3  months  of  the  conference.  This 
conference  signals  the  start  for  a  new  community  of  workers  with  a  common 
Interest  and  this  spreading  of  the  work  Into  other  laboratories  is  the  first 
sign  of  maturity  for  this  new  technology.  A  second  sign  of  maturity  Is  the 
entrance  of  commercial  Instrument  manufacturers  Into  this  field.  We  have 
previously  reported  that  Leltz  Co.  of  Wetzlar,  West  Germany  has  undertaken  a 
program  to  construct  a  prototype  instrument.  We  can  now  report  that  In  Japan, 
Olympus  Optical  Company  has  undertaken  a  strong  effort  In  this  field  and  they 
too  are  following  the  methods  developed  In  our  laboratory.  We  have  received 
copies  of  the  enclosed  bulletin^  as  put  out  by  Hitachi,  also  In  Japan  —  but 
we  know  very  little  of  that  effort  other  than  that  reported  In  this 
publication. 

Our  own  work,  as  reported  below,  consists  of  further  Improvements  In  the 
Instrument.  We  have  extended  the  operation  to  higher  frequencies  (3.6  GHz 
with  a  wavelength  of  0.4  micrometers).  We  have  built  new  hardware  and 
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NEW  TECHNOLOGY 


Noninvasive  Method  for  Visualization 
of  Interior  Region  of  Sample 

» 

Acoustic  Microscope 


In  recent  years,  ultrasonic 
technology  has  been  progressing  by 
leaps  and  bounds.  With  utilization  of 
1  GHz  (Kfi  Hz)  waves  now  a 
possibility,  it  has  become  feasible  to 
develop  acoustic  microscopes  with 
high  resolutions  on  the  sub-micron 
order.  Intensive  research  is  in  fact 
being  directed  toward  the 
development  of  such  microscopes  at 
Stanford  University  in  the  U.S.A.  as 
welt  as  other  research  facilities 
througfiMut  the  world. 

TNs  is  an  area  in  which  Hitachi  has 
put  In  a  lot  of  effort,  too.  In  1978, 
the  company  came  out  with  a 
scanning  traramission  acoustic 
microscope  with  three  micron 
resolution.  Then,  in  the  autumn  of 
1979,  this  same  system  was 
perfected  into  a  reflection  acoustic 
microscope  with  sub-micron 
resolution. 

Fine  Images  Obtained  Using  a 
Micro  Lens  with  High  Numerical 
Aperture 

How  the  reflection  microscope  works: 

1.  A  piezoelectric  film  is  excited  by 
pulse  with  a  high  frequency  carrier 
wave,  causing  it  to  generate  an 
ultrasonic  wave. 


2.  The  ultrasonic  wave  is  sharply 
focused  by  a  spherical  lens  in  the 
sensor  section  and  directed  onto 
the  test  sample  which  is  moved 
two-dimensionaliy  so  as  to  be 
scanned  by  the  ultrasonic  beam. 

3.  The  reflected  wave  from  the 
sample  is  picked  up  by  the  sensor 
section  and  is  then  converted  into 
electrical  signals  for  visual  display 
on  a  cathode  ray  tube. 

Main  Features 

1.  Fine  images  can  be  readily 
obta  .ied  thanks  to  two  new 
developments,  a  lens  system 
capable  of  focusing  the  ultrasonic 
sound  down  to  a  diameter  of  less 
than  one  micron  and  a  scanning 
table  almost  perfectly  free  of 
vertical  fluctuation. 

2.  The  physical  properties  of 
materials  with  complex  structures 
can  be  clearly  determined  since 
the  amplitude  mode  conventionally 
used  in  converting  the  ultrasonic 
wave  4nto  a  visual  image  has  been 
supplemented  by  the  incorporation 
of  a  newly  developed  interference 
mode  which  modulates  the  wave 
into  an  interference  image. 


Cells  Can  Be  Observed  in  Their 
Natural  Living  State 
The  resolution  of  acoustic 
microscopes  has  already  reached 
about  the  same  level  as  that  of 
conventional  optical  microscopes, 
thus  opening  up  numerous  potential 
fields  of  application.  Areas  where 
acoustic  microscopes  are  likely  to 
prove  immediately  useful  include 
medicine  and  biology,  where  they  can 
be  used  to  observe  cells  in  vivo.  They 
are  also  useful  in  physics  and 
engineering,  where  they  offer  a 
means  for  noninvasively  examining 
minute  samples.  One  particularly 
important  application  of  this  type  is  in 
the  inspection  of  the  interior 
condition  of  semiconductor  devices. 
Although  acoustic  microscopes  are 
today  used  almost  exclusively  for 
structure  observation,  it  will,  with  the 
development  of  measurement 
techniques  for  obtaining  quantitative 
physical  parameters,  become 
possible  to  use  them  for  determining 
microscopic  properties  as  well.  This 
will  open  the  way  for  quant.tative 
analysis  of  the  image  characteristics 
of  medical  ultrasonic  diagnostic 
equipment  and  put  the  acoustic 
microscope  on  the  road  toward 
development  as  an  analytical 
instrument. 


Hitachi 


developed  new  software  to  enable  us  to  record  the  Image  tn  digital  form  on 
floppy  discs.  We  have  devised  a  system  for  recording  two  simultaneous  Images 
—  each  with  a  different  frequency.  The  difference  frequency  can  be  varied  as 
we  wish. 

This  we  believe  will  lead  us  Into  the  field  of  acoustic  "spectroscopy", 
a  Important  concept  since  the  evidence  suggest  that  the  objects  that  we 
Investigate  are  sensitive  to  the  acoustic  frequency. 

We  have  done  further  work  on  materials  —  both  inconel  and  manganese 
ferrite. 

We  have  been  advised  by  such  people  as  Marchsesault,  Director  of  Research 
at  Xerox  in  Toronto,  and  Sloan  at  Dupont  that  our  preliminary  work  on  polymers 
Indicates  that  we  are  in  a  position  to  tackle  some  of  the  more  important 
problems  encountered  with  these  materials.  With  this  In  mind,  we  have 
extended  the  theory  for  imaging  of  polymer  samples.*  It  turns  out  that  we 
should  be  able  to  measure  in  a  quantitative  manner  both  the  material  density 
and  the  velocity  of  source  point>by-po1nt  over  the  Image  being  imaged.  This 
will  provide  us  with  a  wealth  of  information  of  use  to  the  polymer  chemists. 

the  work  on  adhesion  of  chrome  on  glass  as  reported  in  the  last 
semiannual  report  has  now  appeared  in  print.  This  paper  has  brought  the 
largest  response  of  any  single  paper  we  have  published.  Thin  films  and 
problems  of  adhesion  seem  to  represent  generic  problems  of  great  interest  to  a 
large  number  of  workers;  Adhesion  is  difficult  to  measure  and  our  method  has 
great  appeal  to  those  working  In  the  field. 

A  strong  effort  has  been  placed  on  photoacoustics  during  the  interval 
covered  by  this  report.  We  now  have  an  argon  laser  properly  modulated  at 
microwave  frequencies  and  Images  are  beginning  to  appear.  The  work  that  we 

*3ee  l^lgs.  3,  4,  and  5  of  Quarterly  Report  marked  G.L.  Report  No.  3173. 
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have  done,  as  reported  here,  concerns  the  theoretical  effort  which  culminates 
In  curves  which  show  the  photoacoustics  response  for  various  combinations  of 
materials.  The  f^otoacoustlc  response  relates  the  mgnltude  of  the  acoustic 
power  generated  per  unit  of  absorbed  optical  power. 

All  In  all,  this  has  been  a  rather  productive  period  v^lch  has  given  us 
new  Inspiration  to  continue  the  momentum  of  our  present  program  of  research. 

II.  NEW  INSTRUMENTAL  ADVANCES 

A  significant  amount  of  effort  was  spent  on  the  design  and  construction 
of  a  new  system  that  allows  digital  control  and  digital  storage  of  the 
Images.  The  system  Is  “stand  alone"  so  that  it  can  be  transported  to 
different  microscopes.  It  incorporates  two  Z-80  microprocessors  —  one  to 
control  the  information  fed  to  the  disc,  and  the  second  to  act  as  buffer. 

This  second  feature  will  permit  us  to  record  the  images  at  a  much  faster  rate 
than  has  been  possible  with  previous  systems.  Our  ability  to  deal  with 
digital  signals  lets  us  work  in  the  realm  that  has  been  called  "digital  image 
processing".  We  can  now  record  a  histogram  and  reprocess  the  signal  by 
adjusting  the  histogram  so  that  it  is  more  or  less  evenly  distributed  over  the 
area  of  the  Image.  It  is  known  that  processing  of  this  sort  can  improve  the 
Image  detail  In  significant  ways.  The  other  system  which  is  now  easy  to 
Implement  Is  the  straightforward  technique  of  differentiating  the  signal  to 
reduce  the  image  contrast  and  bring  out  the  features  such  as  edges.  A 
striking  example  has  been  published  by  Hollis  of  IBM  (see  book^  edited  by 
Ash). 

We  now  have  In  place  a  new  system  that  permits  us  to  record  two 
simultaneous  Images  of  the  same  area.  The  two  Images  can  differ  In  frequency 


by  a  fixed  amount.  With  this  technique  we  make  use  of  two  parallel  channels 
each  with  a  train  of  acoustic  pulses  differing  from  each  other  by  the 
specified  Interval.  The  pulse  trains  are  Interlaced  In  time  and  this  permits 
us  to  separate  them  at  the  output  with  proper  gating  techniques.  We  are  now 
In  the  process  of  extending  the  technique  so  that  we  can  vary  the  frequncy 
differnce  In  a  continous  fashion  and  display  this  In  the  Image.  It  will 
permit  us  to  monitor  and  display  those  features  of  the  object  that  are 
sensitive  to  frequency.  This  flexibility  will  enhance  our  ability  to  extract 
more  and  more  Information  from  the  objects  under  observation. 

III.  ADVANCES  WITH  ACOUSTIC  LENS 

Our  work  has  continued  on  the  theory,  design,  and  fabrication  of  acoustic 
lenses.  This  has  permitted  us  to  build  lens  with  a  smaller  diameter  (now  with 
a  radius  of  16  vm).  This  Is  culminated  In  a  high  Instrument  operating  with  a 
wavelength  of  0.4  micrometers  In  water. 

In  Fig.  1  we  see  an  aluminum  line  Iv  thick  on  a  silicon  substrate  with 
a  2u  sputtered  glass  overlay.  The  acoustic  Images  were  taken  with  lens  of 
radius  equal  to  44  microns  and  operating  at  2.6  GHz.  The  line  Is  15  y  wide. 

It  has  been  Intentionally  damaged  with  a  diamond  Indenter  as  part  of  an 
electromigration  study.  The  optical  micrograph  In  Fig.  1(b)  Is  an  oil 
Ininers Ion  Interference  contrast  picture.  Of  special  Interest  In  the  acoustic 
Image  In  Fig.  1(a)  Is  the  network  of  cracks  In  the  aluminum  line  which  are  not 
evident  In  the  optical  Image  of  Fig.  1(b).  The  dark  area  of  the  Indentation 
In  Fig.  1(a)  Is  seen  to  be  different  In  Its  acoustic  i'esponse  from  the  other 
Indented  areas  and  from  the  rest  of  the  aluminum  line.  A  possible  explanation 
Is  that  the  Indenter  has  broken  loose  the  bond  between  the  aluminum  and  the 


with  quartz  overlay 


silicon  below  or  between  the  aluminum  and  the  quartz  above.  We  know  from  our 
previous  studies  that  adhesion  failure  gives  rise  to  strong  contrast  in 
acoustic  micrographs.  An  alternative  possibility  is  that  the  indenter  has 
removed  the  glass  layer  from  this  area. 

IV.  MATERIAL  STUDIES 

We  have  continued  to  study  the  nickel  alloy.  Inconel.  The  image  of 
Fig.  2  with  its  striking  contrast  is  one  of  the  best  that  we  have  recorded  to 
date.  We  see  there  both  the  grain  boundaries  and  the  curved  lines  and  the 
slip  planes  as  the  straight  lines.  The  slip  planes  are  more  difficult  to 
delineate  with  chemical  etching  procedures,  yet  in  the  acoustic  micrographs 
they  show  up  with  a  contrast  equal  to  that  of  the  grain  boundaries.  The 
important  feature  in  these  images  is  the  strength  of  the  acoustic  reflection 
at  the  grain  boundaries.  The  acoustic  response  in  this  form  tells  us 
something  of  the  elastic  properties  of  the  boundaries.  As  yet  we  have  not 
unraveled  this  source  of  contrast  but  it  continues  to  interest  us. 

A  second  material  that  we  have  studied  during  this  interval  is  a  sample 
of  manganese  ferrite.  It  is  a  material  commonly  used  for  magnetic  recording 
heads  —  a  material  of  some  technical  importance.  Our  comparison  between  the 
optical  and  acoustic  micrographs  shows  that  the  surface  features  appear  in 
both,  but  the  grain  boundaries  appear  only  in  the  acoustic  micrographs.  The 
relative  orientation  of  the  grains  can  be  deduced  from  the  change  in  contrast 
between  adjacent  grains.  The  one  defect  is  seen  as  a  chip  at  the  narrow 
junction  between  three  of  the  grains.  As  a  result  of  this  study,  our  friends 
at  Honeywell  Laboratory  in  Denver  have  provided  us  with  samples  of  actual 
recording  heads.  We  hope  to  be  able  to  present  the  results  of  that  study  in 


INCONEL  ALLOY 


(a)  OPTICAL,  POLISHED 


(c)  ACOUSTIC  Z=  -0.5fi 
2.7  GHz 


(b)  OPTICAL,  ETCHED 


(d)  ACOUSTIC  Z  =  -I.O 
2.7  GHz 


FIGURE  2 


our  next  report.. 


We  mentioned  in  our  introduction  that  the  publication  of  a  paper  on  film 
adhesion  has  attracted  attention.  One  instance  of  this  which  is  of  direct 
interest  to  the  Air  Force  is  the  call  we  received  from  the  Hughes  Laboratory 
at  Santa  Barbara.  They  are  experiencing  problems  with  the  adhesion  of  Lead 
Sulphide  films  to  the  sapphire  substrate.  These  devices  are  used  as  infrared 
detectors.  The  project  is  being  done  for  the  Space  Missle  Division  in  Los 
Angeles.  Dr.  Jim  Knight  at  Santa  Barbara  has  provided  us  with  defective 
detectors  and  we  are  currently  in  the  process  of  characterizing  these  devices 
with  acoustic  microscopy. 


V.  PHOTOACOUSTIC  MICROSCOPY 

From  our  previous  reports  the  reader  will  be  aware  of  our  strong  interest 
in  photoacoustics  as  a  new  method  for  imaging.  In  essence,  it  will  allow  us 
to  monitor  and  image  the  optical  energy  absorbed  by  the  object  that  is  under 
scrutiny.  The  modulated  optical  beam  periodically  heats  the  sample  and 
generates  an  acoustic  wave  in  the  liquid  surrounding  the  sample.  This  wave  is 
picked  up  in  the  acoustic  lens  and  the  image  is  formed  with  mechanical 
scanning  as  in  the  conventional  acoustic  microscope.  The  new  system  is 
operating  with  acoustic  waves  near  900  MHz.  We  are  now  recording  our  first 
traces  with  a  good  signal-to-noise  ratio. 


Our  plan  is  to  present,  in  this  section,  the  results  of  the  theoretical 
calculations  that  we  have  been  performing  in  order  to  predict  that  magnitude 
of  the  acoustic  wave  that  can  be  expected  for  a  given  amount  of  absorbed 
optical  power. 

1.  Scanning  Photoacoustic  Microscopy 

In  the  photoacoustic  effect  absorbed  light  is  converted  into  sound  by 

the  thermal  expansion  of  the  absorber  and  its  surroundings.  It  is  the  basis 

of  photoacoustic  spectroscopy,  where  the  received  acoustic  signal  detected 

by  a  microphone  is  measured  as  a  function  of  the  optical  wavelength  of  light 

3  4 

absorbed  by  a  solid,  liquid,  or  gas.  ’  It  can  be  used  to  measure  very  small 

5  6 

concentrations  of  materials,  or  very  small  absorption  coefficients.  *  Early 

24 

flaw  detection  results  by  von  Gutfeld  pointed  to  the  possibility  of  micros¬ 
copy  based  on  photoacoustics. 

In  the  work  to  be  described  here  we  will  first  review  the  work  that  was 
done  prior  to  this  reporting  interval  in  order  to  remind  the  reader  what  was 
done  and  to  prepare  for  the  next  report  where  we  expect  images  of  superior 
quality  as  recorded  on  the  new  machine.  In  the  later  section  on  theory  we 
report  the  work  that  was  done  during  this  interval.  We  see  there  that  the 
results  are  presented  in  forms  of  curves.  These  are  most  useful  in  pre¬ 
dicting  the  acoustic  output.  Figure  3  shows  the  experimental  setup  that 
was  used  initially.  It  has  now  been  replaced  with  an  argon  laser  that  is 
modulated  with  an  acousto-optic  modulator  at  a  single  frequency.  An 
acoustic  lens  of  radius  of  200  pm  was  used  as  the  receiver. 
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The  confocal  arrangement  means  that  we  collect  a  substantial  fraction 
of  the  generated  acoustic  energy.  This  Is  essential  to  overcome  acoustic 
losses  In  the  water  used  as  coupling  fluid  from  sample  to  lens.  The  use  of 
the  acoustic  lens  ensures  that  our  output  acoustic  signal  Is  collected  only 
from  the  heated  spot  on  the  sample.  A  potential  advantage  of  the  high 
acoustic  frequency  Is  the  ability  to  generate  acoustic  efficiently  with 
short  modulated  optical  pulses.  We  shall  see  later  that  the  strength  of  the 
photoacoustic  signal  Is  a  strongly  Increasing  function  of  the  optical  modu> 
latlon  frequency  for  some  materials. 

The  sample  was  scanned  with  a  scanner  from  a  conventional  acoustic 
microscope,  modified  to  scan  at  1  Hz  fast  scan  rate.  Position  was  detected 
In  both  X  and  y  directions  perpendicular  to  the  acoustic  lens  axis. 

These  signals  controlled  the  position  of  the  CRT  spot.  The  Intensity  of  the 
cathode  ray  display  was  proportional  to  the  amplified  and  detected  acoustic 
power  In  a  pulse. 

The  photoacoustic  signal  Is  generated  by  the  intensity  modulation  of 
the  absorbed  laser  light.  Later  we  will  discuss  the  theory  of  photoacoustic 
generation  and  develop  a  detailed  model  predicting  photoacoustic  output  for 
given  laser  input  for  different  geometries.  For  now  we  note  that  the  acoustic 
wave  equation  In  the  presence  of  thermal  expansion  for  a  semi-infinite  Iso¬ 
tropic  solid  1s:^ 


OBJECT 


FIG.  3— The  photoacoustic  setup. 


where 


T  ■  stress 

*  sound  velocity 
B  «  bulk  modulus 

B  •  linear  coefficient  of  thermal  expansion 
Cji  *  stiffness 
e  ®  temperature  rise. 

By  inspection  of  this  equation  we  can  see  the  right  hand  side  is  the  driving 
or  source  term  for  the  stress  wave  T(z,t)  .  One  can  also  see  that  the  fre¬ 
quency  spectrum  of  the  temperature  variation  e(z,t)  will  determine  the 
spectrum  of  the  generated  acoustic  waves. 

Now  the  temperature  rise  6  is  proportional  to  optical  power  absorbed, 
providing  the  absorption  coefficient  o  is  not  a  function  of  temperature. 
Thus  the  acoustic  spectrum  is  determined  by  the  spectrum  of  the  modulation 
of  optical  power  absorbed  by  the  sample. 

We  also  notice  that  in  terms  of  power,  with  9  «  ^optical 

''.coustic  *  ’’acoustic  ”  '’’optical  ' 

In  order  to  produce  acoustic  images  with  this  system  the  system  was 

confocally  aligned,  so  that  the  optical  lens  and  the  acoustic  lens  are 

focused  on  overlapping  spots.  This  is  a  disadvantage  of  the  present  system 

as  compared  to  other  systems,  which  typically  use  a  microphone  or  hydrophone 

in  a  gas  or  liquid  cell  as  the  acoustic  receiver.  There  the  difficult  and 

precise  alignment  of  the  system  is  not  required. 

The  resolution  we  expect  in  this  confocal  imaging  system  is  better  than 

we  get  from  one  lens  alone.  Basically,  we  expect  the  overall  imaging  system 

transfer  function  to  be  the  convolution  of  the  transfer  functions  of  each 
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lens.  This  means  that  higher  spatial  frequencies  in  the  object  are  passed 
by  the  imaging  system  and  by  the  criteria  discussed  elsewhere.  This  means 
greater  resolution.  A  complicating  factor  here  is  that  acoustic  waves  will 
be  produced  from  a  source  larger  than  the  optical  spot  size  due  to  thermal 
diffusion.  Thus  at  840  MHz  the  thermal  diffusion  length  in  gold  is  0.23  p. 
If  we  have  an  optical  spot  size  of  2.0  p  we  will  have  a  thermal  spot  size 
(where  the  temperature  is  modulated  at  840  MHz)  of  ~  2.5  p.  The  square-law 
nature  of  the  photoacoustic  generation  process  acts  to  offset  this  by 
shrinking  the  effective  spot  size. 

To  produce  confocal  alignment  of  optical  and  acoustic  lenses »  the 
acoustic  lens  acted  as  the  lens  of  a  reflection  acoustic  microscope  to 
image  a  scanned  finder  grid.  A  scanning  optical  reflection  micro¬ 
scope  was  constructed  by  picking  off  with  a  beam  splitter  a  portion  of 
the  light  returning  through  the  optical  lens  from  the  other  side  of  the 
same  finder  grid.  The  optical  reflection  microscope  image  was  formed 
by  focussing  this  returning  light  onto  a  photodiode  whose  output  con¬ 
trolled  the  intensity  of  the  display.  The  same  X  and  Y  position 
signals  from  the  scanner  were  used  to  position  the  spot  on  the  face  of 
the  display.  By  adjusting  the  positions  of  the  optical  and  acoustic 
lenses  with  micrometers  on  their  respective  stages,  one  could  bring 
the  two  reflection  images  -  optical  and  acoustic  —  into  approximate 
coincidence.  One  could  improve  alignment  by  forming  a  difference 
image  by  subtracting  acoustic  and  optical  signals  in  a  differential 
amplifier.  Finally  input  laser  power  was  increased,  the  acoustic  system 
was  converted  to  receive-only  mode,  and  one  looked  for  photoacoustic 
pulses  on  an  oscilloscope  triggered  by  the  laser  shot.  The  signals  were 
typically  not  visible  without  further  search  by  systematically  'walking' 
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one  set  of  stage  micrometers  around  the  position  of  rough  coincidence 
already  arrived  at.  Once  located,  a  photoacoustic  signal  could  be  kept 
In  sight  for  several  hours.  The  Images  could  be  made  over  various  areas 
of  the  sample. 

To  produce  the  Image,  the  detected,  sampled  and  held  photoacoustic 
pulse  supplied  the  Intensity  control  for  the  display,  while  sample  posi¬ 
tion  sensors  X  and  Y  controlled  CRT  spot  position. 

The  earliest  Images  produced  with  this  technique  were  anything  but 
nondestructive.  In  particular.  Images  of  a  1  mil  (25  y)  thick  A1  foil 
and  a  Cr/AI  finder  grid  pattern  were  imaged  with  this  technique  and  these 

O 

images  have  been  published  elsewhere.  Damage  to  the  samples  due  to 
the  high  temperatures  reached  during  imaging  was  severe.  Holes  were 
burned  through  the  1  mil  A1  foil,  and  both  A1  and  Cr  in  the  grid  pattern 
were  severely  damaged.  The  melting  temperature  of  Aluminum  Is  660  C; 
that  of  Chromium  Is  1857  C.  We  conclude  that  local  temperature  changes 
of  2000  C  or  more  were  produced  In  this  first  imaging  effort. 

In  a  second  round  of  imaging  with  this  photoacoustic  system  the 
scanning  and  positioning  system  was  improved.  In  Fig.  4  we  show 
optical  and  photoacoustic  images  of  a  gold  finder  grid  1500  A  thick 
evaporated  on  a  microscope  cover  slip.  1.8  y  of  sputtered  glass  (X/4 
at  840  MHz)  was  deposited  on  top  of  the  gold,  on  half  of  the  finder 
grid  pattern,  including  the  letter  *V*  as  shown.  The  sputtered  glass 
layer  was  Intended  to  serve  as  a  quarterwave  acoustic  transformer  from 
gold  Into  water.  It  was  observed  to  function  as  such,  since  the  portion 
of  the  finder  grid  covered  with  the  sputtered  glass  was  observed  to 
have  dB  higher  photoacoustic  signal  than  the  portion  of  the  grid  not 
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X  =  1.83/i. 

FIG.  4  — Photoacoustic  image  of  chrome  on  glass. 
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covered  with  this  glass. 

In  the  optical  reflection  microscope  picture  of  the  letter  'V’, 
Individual  glass  particles  from  the  sputtering  process  are  visible  on 
the  gold  pattern.  Many  of  these  are  recognizable  as  dark  areas  In  the 
photoacoustic  Image  as  well.  Presumably  they  act  as  scattering  centers 
for  the  photoacoustic  pulses. 

At  the  bottom  of  the  circular  band  surrounding  the  letter  'V  Is 
an  Imperfection  in  the  pattern  visible  In  the  optical  micrograph  which 
Is  not  visible  in  the  photoacoustic  Image.  It  Is  a  small  area  of  metal 
which  was  evaporated  onto  photoresist  In  the  process  of  laying  down  the 
grid  pattern,  it  did  not  break  away  cleanly  in  the  liftoff  pro¬ 
cedure.  This  area  of  metal  is  still  attached  to  the  well  adhered  grid 
pattern,  but  is  not  itself  well  adhered  to  the  glass  substrate.  The 
area  in  question  was  evaporated  onto  photoresist,  which  was  dissolved 
away  during  liftoff.  It  Is  laying  on  the  glass  substrate  and  is  covered 
with  the  sputtered  glass. 

The  lack  of  adhesion  is  shown  by  the  photoacoustic  picture,  which 
simply  does  not  show  this  area  at  all.  The  photoacoustic  picture  was 
taken  before  the  opcical  picture.  If  we  assume  this  piece  of  metal 
was  left  on  at  liftoff,  it  must  have  been  physically  present  in  the 
sample  when  the  photoacoustic  image  was  made.  .  It  must  have  been 
the  different  boundary  conditions  at  the  gold-glass  Interface  which 
resulted  In  no  photoacoustic  signal.  As  we  shall  see,  the  theory  of 
photoacoustic  generation  predicts  that  light  absorbed  at  a  free  sur¬ 
face  Is  much  less  efficient  In  generating  a  photoacoustic  signal  than 
a  constrained  surface.  The  non-adhering  metal  area  we  have  been  dis¬ 
cussing  Is  seen  to  act  as  a  free  surface. 


The  optical  micrograph  in  Fig.  4  was  taken  after  the  photoacoustic 
Imaging.  There  is  no  visible  damage  in  the  optical  picture.  In  other 
pictures  taken  of  other  areas  of  the  same  glass/gold/glass  grid  it  was 


found  that  with  optical  power  approximately  twice  that  in  Fig.  4,  the 
gold  started  to  melt.  The  melting  point  of  gold  is  1064  C.  We  con¬ 
clude  that  temperature  changes  of  several  hundred  degrees  C  were 
encountered  in  the  photoacoustic  imaging  of  Fig.  4. 

In  Figs.  5  and  6  we  see  photoacoustic  images  of  an  integrated 
circuit.  This  circuit,  a  test  silicon  on  sapphire  device  manufactured 
by  RCA,  was  illuminated  by  the  laser  pulses  through  the  sapphire  sub¬ 
strate.  The  laser  light  was  absorbed  at  the  bottom  of  the  layers 
making  up  the  integrated  circuit  structure.  Sound  propagated  from 
the  front  (top)  side  of  the  devices  to  the  acoustic  lens.  The  photo¬ 
acoustic  signal  is  absent  where  there  are  no  optical  absorbers  present. 
Different  image  intensity  is  seen  for  different  types  of  device  layers. 
A  feature  worthy  of  note  is  that  areas  of  uncovered  silicon  (visible 
as  lightly  covered  areas  in  the  optical  pictures)  give  much  weaker 
photoacoustic  response  than  areas  of  a  aluminum  or  aluminum-covered 
silicon  (black  in  the  optical  pictures). 

2.  Theory  of  Photoacoustic  Generation 

To  predict  the  photoacoustic  response  we  will  adopt  a  plane-wave 
model.  The  incoming  in tensity -modulated  light  will  be  absorbed  and 
raise  the  temperature  of  the  absorbing  region.  The  heat  diffusion 
process  will  determine  temperature  rise  as  a  function  of  position. 

This  temperature  function  is  used  as  the  driving  term  in  the  acoustic 
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field  equations,  which  are  solved  to  find  the  acoustic  power  in  the  far 

9 

field.  The  method  of  coupled  modes  is  used  to  solve  for  the  generation 
from  the  distributed  source  region. 

Figures  7  and  8  show  the  specific  models  we  are  postulating.  In 
Fig.  7  light  enters  as  a  plane  wave  of  intensity  through  trans¬ 
parent  medium  I  .  It  is  absorbed  in  medium  II,  the  light  intensity 
varying  as  1^6”°*^  .  Heating  and  resultant  thermal  expansion  of  the 
heated  boundary  region  generates  acoustic  waves  that  we  detect  in  the 
far  field  of  region  I  or  region  II.  We  eventually  assume  the  incoming 
light  is  intensity  modulated  with  angular  frequency  w  ,  and  the  re¬ 
sultant  acoustics  are  at  the  same  frequency.  In  Fig.  8  we  make  region  II 
a  layer  of  thickness  L  and  there  the  light  is  attenuated  as  • 

Region  III  is  transparent.  We  will  solve  for  the  temperature  distribu¬ 
tion  as  a  function  of  position,  and  use  this  as  the  driving  term  in  the 
acoustic  field  equations.  We  will  solve  these  with  coupled  mode  vari¬ 
ables  for  the  far-field  acoustic  strain  and  power. 

3.  Two-Layer  Model;  Heat  Flow 

The  theory  of  heat  conduction^®  tells  us  that  heat  will  flow  when¬ 
ever  we  have  a  temperature  gradient  in  space.  For  the  case  of  one¬ 
dimensional  heat  flow,  we  will  have: 

f  =  -k||  (5.1) 

2 

where  f  is  the  heat  flux  in  watt/m  ,  0  is  the  temperature,  and  K 


REGION  I 


a|.|(2) 


\ 


REGION  IL 


!  ''•si  = 


z  =  0  “Utz) 


(o) 


✓ 

I  \ 


s 


fl.= 


0»  = 


z  =  0 


(b) 


FIG.  7.  (a)  Production  of  the  thermal  wave,  two-region  case, 

(b)  Temperature  profile  for  the  two-region  case. 
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FIG.  8v  (a)  Thermal  waves  in  the  three-region  case, 
(b)  Temperature  profile  three  regions. 


is  the  thermal  conductivity  in  (watt/km).  If  the  heat  flux  changes 
with  position  z  ,  there  is  heat  being  deposited  or  generated  in  the 
medium.  The  rate  of  heat  absorption  in  a  volume  element  AzAA  is: 

-(|£az)aa 

where  the  negative  sign  accounts  for  the  fact  that  if  3f/3z  is  a 
negative  number,  net  energy  is  being  absorbed  in  the  medium,  if  we 
assume  a  flow  of  heat  toward  positive  values  of  z  . 

The  medium  may  also  be  heated  by  absorption  of  optical  energy. 

For  light  traveling  towards  +z  where  the  intensity  varies  as 

'  'o® 

we  will  have  energy  absorbed  by  the  medium  equal  to: 

-(|Jaz)aa 

for  a  volume  element  AzAA  .  Both  the  change  in  heat  flux  and  the  ab- 
soprtion  of  optical  energy  will  cause  a  change  in  the  temperature.  For 
a  volume  AzAA  this  will  correspond  to  rate  of  energy  absorption 

AC  II  AZAA 

where  p  =  density;  C  =  specific  heat  .  Equating  heat  absorbed  to  heat 
required  for  temperature  change  0  ,  we  find  that: 


We  now  assume  sinusoidal  excitation*  and  this  gives: 


•  -jwce  ♦ 


where  (5.4)  simply  restates  the  definition  of  heat  flux  f 


(5.3) 


(5.4) 


The  standard  method  of  solving  these  equations  for  the  temperature 

♦ 

0(2)  Is  to  eliminate  the  variable  f  .  One  is  left  with  the  diffusion 
equation,  and  there  Is  a  large  literature  devoted  to  solving  this  equa¬ 
tion  under  a  variety  of  boundary  conditions. However,  to  point  out 
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the  advantages  of  coupled  modes,  we  will  use  this  method  to  solve  for 
the  temperature  distribution. 

We  multiply  Eqs.  (5.4)  by  and  add  it  to  Eq.  (5.3): 


2*  /  jwpCK  \ 

■ii'-rj 


+  al  e"*^ 
0 


(5.5) 


Now  let: 


JoipCK 


-  (1  +  j) 


(5.6) 


and  we  have: 


—  (f  +  Z*0) 


-  (f  +  Z^.6)  +  al  e‘ 

^  t  0 


(5.7) 


or,  letting 


(5.8) 


we  have 


.  ^  ai  +  al  e-«" 


(5.9) 


Similarly,  if  we  multiply  Eq.  (5.4)  by  and  subtract  from 
Eq.  (5.3),  we  get 


;,t  .  T  -az 
-  *  —  a  +  aKe 

az  K  ‘  ° 


(5.10) 


where 


f  -  Z^e  . 


(5.11) 


The  advantage  of  transforming  to  normal  mode  variables  is  that 
we  have  uncoupled  the  equations  without  having  to  go  to  a  second  order 
equation.  Equations  (5.9)  and  (5.10)  can  now  be  manipulated  to  a  form 
that  is  directly  integrable.  We  let: 

a*  =  d^  e"^^ 


(5.12) 


where 


t  -  •  (1  +  J) 


(5.13) 
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Then  Eqs.  (5.9),  (5.10)  become: 


^  .  .r 

(5.14) 

^  -  ar 
>2  » 

(5.15) 

We  now  apply  these  equations  to  the  two  regions  I  and  II.  In 
region  II  we  will  have,  by  simple  Integration  of  (5.14): 

dJgCz)  -  d+2(0)  =  ol^  /  e  ^  d; 

''o 

where  c  the  diminy  position  variable  of  Integration.  Integration 
gives: 


d;2(2)  - - Je  ?  -l)  +  d+2(0)  (5.16) 

01  -  T,  \  / 


Similarly,  for  the  negative-traveling  mode  In  region  II, 

.it  /  X  .t  /  ,  /  -(a+T-); 

d_2(«)  -  d_2(z)  ■  al^y  e  ^  dc  (5.17) 

z 

Since  we  have  no  heat  flow  from  +»  ,  d^2(”)  “  0  »  *nd 

d^2(z)  -  - 


Ol^  -(o+T,)z 
e 


X.e  ‘2 

O  +  To 


(5.18) 


We  can  now  put  these  back  In  terms  of  aj  and  to  get: 


t  .  "«Z\  -ToZ 

aj  (z)  -  + - a-  (e  2  .  e  )  +  a  2(0)  e  2 

O  -  To  '  / 


(5.19) 
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(5.20) 


F 


a!, (2) 


al 


o_  g-oz 


O  ♦  T, 


In  region  I  we  have  no  sources  of  heat,  and  hence  Eqs.  (5.10)  and 
(5.11)  simplify  to: 


and 


da 


_+l 

dz 


®+l 

•^1 


t 

-1  • 


Since  a^j(-“)  *  0  ,  we  will  have 


a+^(z)  =  0 
and 

a^i(z)  -  a^j(0)e^^^ 


(5.21) 


What  remains  Is  to  match  boundary  conditions  at  the  interface 
Z  «  0  .  There  we  must  have  continuity  of  the  temperature  0  and  the 
heat  flux  f  .  For  the  latter  (f)  we  must  have; 

«n(0)  +  a^j(O)  =  aJgCO)  +  a^2(0)  (5.22) 

For  the  former  (8)  we  will  have; 

fi2i(ali(0)  -  a!i(0))  “ 

Substitution  of  our  results  for  ajj(z)  ,  a^j(z)  ,  aJgCz)  and  a^2(^) 
Into  these  gives: 

-  27  - 


(5.24) 


»!l(0) 


(691  +  l)(a  +  T,) 


•^(o) 


(1  » 

(1  *  421^“  +  ^2) 


(5.25) 


where  we  have  let  21  *  * 

Using  Eqs.  (5.19),  (5.20) «  and  (5.21)  we  can  write  for  the  tem¬ 


perature  d1 strl buti ons : 


61(2)  =  0ije 


89(2)  =  090®  +  099® 


(5.26) 


(5.27) 


where: 


Q  ,  _ 2 _ 

Zj(l  +  62i)(a  +  T2) 


oIq  (1  +  62i(a/T2)) 


(5.28) 


(5.29) 


If  we  let 


Wi 


(5.30) 


(5.31) 
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It  nuty  be  shown  that: 


(5.32) 


®22  '  ”  ®02 


(5.33) 


Again  we  have  defined: 


^f.2  * 

(1  +  j) 

V  2  ^1.2  ‘'l.Z  '^1.2 

(5.34) 

^1.2  • 

n.2 

'^1,2 

(5.35) 

‘21  ■ 

/P2  C2  K2‘ 

^1  ^1  ^1 

(5.36) 

To  calculate  the  maximum  temperature  excursion  of  the  sample,  we 
note  that  the  peak  temperature  will  be  at  z  =  0  ,  where 


©iCO)  *=  Ojj  *  $2(0)  =  OgQ  +  022  • 

In  Table  5.1  are  listed  some  calculated  peak  ac  temperature  excursions 

at  the  Interface  z  =  0  for  a  variety  of  combinations  of  materials  for 
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a  100  MHz  modulation  frequency,  with  10  W/m  »  10  mU/u  absorbed.  We 
note  the  remarkable  similarity  In  results  for  various  types  of  medium  1 
used  with  the  metals  and  silicon.  This  Is  because  the  metals  are  better 
heat  conductors  than  the  first  media  (pyrex,  sapphire,  water),  so  the 
first  medium  Is  of  minor  consequence. 


TABLE  5.1. 

Maximum  Temperature  Excursion  at  100  MHz  for  Optical 

Intensity  Absorbed 

I  -  10^ W 

»  10  mW//. 

Medium  1 

Medium  2 

u 

-l<*) 

a,  m 

Temperature  0^^ 

Pyrex 

At 

1.12  E8 

19  C 

Sapphire 

At 

1.12  E8 

16  C 

Water 

AA 

1.12  E8 

19  C 

Pyrex 

Au 

7.60  E7 

18  C 

Sapphire 

SI 

3.29  E6 

16  C 

Pyrex 

Hater 

1.00  E5 

0.37  C 

Pyrex 

Water 

1.00  E6 

3.7  C 

4.  Two-Layer  Model;  Acoustic  Generation 

Having  solved  for  the  temperature  profile  excited  by  light  inten¬ 
sity  modulated  at  radian  frequency  m  ,  we  can  proceed  to  calculate  the 
generation  of  acoustic  waves.  The  coupling  from  a  fluctuating  tempera¬ 
ture  distribution  to  propagating  sound  wave  is  via  the  coefficient  of 
thermal  expansion. 

For  a  one-dimensional  isotropic  model  ,  Hooke's  law  takes  the 

form: 

T  •  CjjS  -  Be'0  (5.37) 

where  T  »  stress  ,  S  “  strain  »  Cjj  “  stiffness  ,  B  »  bulk  modulus  , 

P'  ■  linear  thermal  expansion  coefficient  ,  0  =  temperature  change. 
Bulk  modulus  B  rather  than  c^^  Is  used  In  the  second  term  because 
by  the  nature  of  thermal  expansion  the  stress  It  causes  Is  a  hydro- 

'  static  pressure  Increase  rather  than  a  directed  traction  force. 


For  sound  propagation  in  one  dimension  (z)  in  an  isotropic  medium 


we  have: 


(S.38) 


where  u  *  particle  velocity.  This  is  just  the  time  derivative  of 
the  definition  of  strain.  Substituting  in  (5.37)  we  have: 


M  -  JLii  ft  li 
az  "  Cjj  at  ‘  3t 


where  e  =  Bg'/Cjj  .  Newton's  second  law  for  this  case  is: 


(5.39) 


aT  _  du 

Si- 


(5.40) 


Equations  (5.39)  and  (5.40)  are  two  coupled  equations  in  the  variables 
T  (stress)  and  u  (particle  velocity)  .  Letting  «  pV^  ,  where 
Z  ■  acoustic  impedance,  V  =  sound  velocity,  we  define 


T  +  Z  u 
0 


a  ■  T  -  Z^u 
-  0 


(5.41) 

(5.42) 


where  a^  and  a_  (with  no  superscript  as  in  a^,  a^)  are  positive 
and  negative  traveling  acoustic  normal  modes. 

We  assume  that  temperature  waves  and  sound  waves  are  time  harmonic 


and  hence  a/ at  jm  .  Then 


It  •  -j«pu 


liL  ■  .  J<u 


T  “  ju06 


(5.43) 


(5.44) 
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To  transform  to  normal  modes  we  multiply  Eq.  (5.44)  by  and  add  to 
Eq.  (5.43): 


^  (T  +  Z^u)  »  -j  ^  T  -  j<apU  -  .j«e  z^e 
^11 

With  Cjj  ■  pVg  and  k  *  m/V^  : 

^  (T  +  Z^u)  -  -jk(T  +  Z^u)  -  jkCjjPe 
or: 

Similarly,  if  we  multiply  Eq.  (5.44)  by  Z^,  and  subtract 
Eq.  (5.43)  we  get: 

^(T-Z„u)  .  (T-Z„u) 

s 

or 

We  make  a  change  of  variables  to  carry  out  Integration: 

and  a.  - 

This  changes  Eqs.  (5.47)  and  (5.49)  to: 


(5.45) 


(5.46) 


(5.47) 

It  from 


(5.48) 


(5.49) 


and 


—  •  +Jkc,,pee  (5.51) 
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If  we  now  assume  6  Is  different  from  0  over  a  region  extending 
from  Zj  to  Z2  we  can  Integrate  these  equations: 


d^(z)  -  d^(2j)  «  -jkCjjB  f  ee'^'^'^^  dc 

(5.52) 

•r 

21 

^2 

d_(z2)  -  d_(z)  *  +jkCjjB  J'  6e"^*^^  dt 

(5.53) 

which  yield: 

a+(z)  =  -jkCjjBe^  J"  d;  +  a^(Zj)e  ^  ^  (5.54) 

^1 

and 

a,(z)  *  -jkCjjPe'^^’^^  J  d;  +  a_(z2)e  (5.55) 

z 

Let  us  now  return  to  the  specific  case  of  the  two  semi -Infinite 
media  meeting  In  the  plane  z  =  0  ,  with  optical  energy  Incident  from 
the  left  In  transparent  medium  I  and  absorbed  in  medium  II.  In  medium 
II  we  have  previously  found  the  temperature  in  Eq.  (5.27)  to  be: 


89(2) 


-T«Z  -OZ 

®20®  *  ®22® 
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the  boundary  conditions  can  be  written  as: 


+  a_j(0)  =  (5.60) 

4jj(a^j(0)  -  a.j(O))  •  a^,2(0)  -  a.j(O)  (5.61) 

Where  we  have  defined  Agj  =  202^^01  *  Substituting  the  values  for  the 
various  a's  from  Eqs.  (5.56),  (5.57),  (5.58),  and  (5.59),  into  Eqs. 
(5.60)  and  (5.61),  we  can  solve  these  equations  for  the  remaining  previous 
unknowns,  results  are: 


2Z 


02 


^0l‘^^02 


1^1+  jk^  j 


p0l"^02\  j 

'  W^02' 


jk2Cii2B2 


’22 


’20 


a  +  jk2  ^2  jk2 


(5.62) 


a.i(O)  =  - 


2Z 


01 


^0l'*^^02 


jk2Cii2B2  +  —hl-\ 

^  ^  \«  +  jk«  Tp  +  jkJ 


/^or^o2\ 

»  • 

^'^l^l'^lll®!! 

'^0l'*^^02' 

(5.63) 


Finally,  in  the  far-field  region  where  the  sound  will  be  observed, 
Eqs.  (5.56)  and  (5.57)  simplify  somewhat  because  in  the  far  field  all 
terms  containing  exponentials  with  real  parts  go  to  0,  leaving  only 
terms  containing  .  Then  in  this  limit  we  will  have: 


-  35  - 


^  ^  ^  jk2  -  a  jk2  -  T2  \ 


Jk«  c«  1 4  4 

a  (z)  e  1.11  +  a 

jkj  -  Tj 


.l(0)Je  ‘ 


(5.65) 


If  we  substitute  in  for  a^j(O),  a_j(0)  f  and  get  e^j  ,  we  get,  after 


some  rearrangement. 


= 


2jki3iCiiiA2i  I  -  ij  621  •i'^2^112^2 

Tj(1  +  A2i)(1  +  62i)  (1  +j  — )  T2 


1-j-^  Vl  +  Agj/Vl  +  j 


1  +  60,  — 


1  +  6, 


l-j—  \l  +  A2j/  \l  +  0 


-OkgZ 


’^•'1^111^1^21  (^  '  ^)  ’  ^2l\  /  ^  \  1  1 

Tid  t  «jj)  \i  ♦  17775/  7775] 


(5.66) 


2jk2Cii262  (1  +  621  f^)  T2/a 

V' * ''21)  [  (1  ♦  «2j)  (i775j  7775 


(5.67) 


These  expressions  are  the  complete  solutions  for  the  far-field 
acoustic  modes,  and  are  rather  complicated  as  they  stand.  We  would 
like  to  consider  them  in  limiting  cases  of  interest.  First,  suppose 
we  have  strong  optical  absorbtion  (a  »  Tg  ;  a  »  kg)  and  moderate 
diffusion  (x^  »  k.  ;  i=l,2).  We  will  keep  terms  of  the  form 
AgjCkg/xg)  ,  because  may  be  large.  Then  substitution  of  these 
limits  and  substantial  rearrangement  gives: 


-jk.2 

a_j(2)e  ^ 


®+2 


jkgZ 


^^0^21 


$2^2] 


(1  +  AgiKl  +  6gj)  L  Cg  J 


1  + 


{'■•it)  (a  (If 


2'o«21 


^s2^2 

(l  +  Agj)(l  +  6gj)  [  Cg 


1  + 


'-a  0(3'" 


(5.68) 


(S.69) 


where  ic^.  »  thermal  diffusivity  =  {K^./p^c^)  . 

We  can  simplify  even  further  by  assuming  that  medium  1  has 

properties  identical  to  medium  2  (except  that  1  is  transparent 
while  2  absorbs  light).  Then  we  find  that: 


jk.,z  -jk.z 

=  a^j(z)e 


-  y  [2  +  j(k/T)J  ,  (5.70) 
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In  the  low-frequency  limit  as  k/r  0  , 


jk-2  -jk.z 

a^2(z)e  *  a.l(z)e 


(5.71) 


This  reduces  the  photoacoustic  conversion  problem  to  its  simplest 
form.  There,  are  no  reflections  of  acoustic  or  thermal  waves. 

We  see  that  the  simple  expression  [V^P/C]  gives  the  conversion 
of  light  intensity  to  acoustic  stress, wave  amplitude  directly.  In 
terms  of  acoustic  intensity,  we  have  in  this  limiting  ideal  case: 


ac 


1 

K] 

■p- 

8 

.  P  . 

C 

^0  =  ’^^o 


(5.72) 


2 

The  proportionality  constant  tj  between  and  is  seen 

as  an  intrinsic  conversion  factor  characterizing  a  material.  It 
is  the  efficiency  of  conversion  of  light  to  sound  in  the  absence  of 
effects  due  to  mismatch  in  the  properties  of  materials. 

In  Table  5.2  we  list  the  thermal  and  acoustic  properties  of  se¬ 
lected  materials.  The  intrinsic  conversion  factor  n  for  these  nateri 
als  is  listed  in  the  table.  We  note  that  the  only  gas  in  the  table, 

STP  air,  is  by  far  the  most  efficient  converter,  owing  to  its  high 
expansion  coefficient  and  low  density.  The  liquids  in  the  table  are 

more  efficient  than  the  solids,  again  owing  to  larger  expansion  coef¬ 
ficients.  The  conducting  solids  listed  are  more  efficient  than  the 
semiconductors  and  insulators,  although  one  could  find  metallic  sub¬ 
stances  such  as  Invar,  which  would  have  conversion  efficiency  zero 
owing  to  a  vanishing  expansion  coefficient. 
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Returning  now  to  the  case  of  two  different  media  In  the  limiting 
case  of  «  a.  ^2  «  o,  as  expressed  In  Eqs.  (5.68), 

(5.69)  we  note  that  the  positive  and  negative  traveling  modes  a_^ 
and  a^2  ®'’‘*  Independent  of  frequency  (provided  the  fre¬ 
quency  Is  sufficiently  small  enough  that  ^ 

«  1).  The  acoustic  intensities  carried  away  to  the  far 
fields  by  these  modes  are  unequal,  however,  since 


The  acoustic  intensity  will  be  inversely  proportional  to  the  acoustic 
'Impedance.  Me  are  postulating  conditions  here  corresponding  to  a  thin 
plane  transducer  at  the  boundary  between  two  half-spaces.  The  ratio  of 
the  acoustic  intensities  radiated  into  the  two  regions  in  this  case  is 
just  what  we  have  found. 

At  frequencies  where  A2j(k2/T2)  or  (kj/A2jTj)  cannot  be  neg¬ 
lected  compared  to  1,  the  acoustic  mode  amplitudes  In  Eqs.  (5.68) 
and  (5.69)  become  frequency-dependent.  Assume  for  example,  that 
A2j(k2/T2)  ^  1  ,  but  (kj/A2jTj,)  «  1  .  This  obviously  may  happen  if 
&22  »  1  .  As  an  example,  if  medium  2  ®  gold,  medium  1  =  water, 
and  A22  *  39.1  f  then  Eqs.  (5.68)  and  (5.69)  become: 


a,l(2)e 


(l+A2i)(l  +  «2p 


^2 


"  40  « 


X 


(5.74) 


jk,z  2i«-, 

a+2(2)e  ^ - — - 

2  (1  +  A2i)(1  +  «2i)  I  ^2  J 

■[■•(S)t7r-(^)(‘4V] 

In  this  case  the  negative- traveling  mode  is  seen  to  remain  flat  as  a 

function  of  frequency,  whereas  the  positive-traveling  mode  has  the 

form  A  +  (l  +  j)Ba>  '  where  A  and  B  are  real  expressions.  The 

2  1/2  2 

acoustic  intensity  then  has  the  form  I  _  «  A  +  2ABto  '  +  2B  w  . 

cC 

In  Fig.  9  are  plotted  the  intensities,  Ij  ,  and  I2  of  acoustics 

generated  in  media  1  and  2  for  the  combinations  methanol -aluminum 

(^21  -  21)  ;  water-aluminum  (a2j  =  11)  ;  pyrex-aluminum  (Agj  =  1.33)  ; 

sapphire-aluminum  {A2J  =  .39)  .  Aluminum  is  medium  2  in  each  case.  The 

full  expressions  [Eqs.  (5.66)  and  (5.67)]  with  no  approximation  are  used. 

A  listing  of  the  program  used  is  given  in  Appendix  C.  The  ordinates  are 

2 

acoustic  far-field  intensity  in  watt/m  for  optical  absorption  of 
2 

1  watt/m  .  The  main  features  of  those  graphs  are  in  accord  with  our 
previous  discussion.  Ij  is,  in  each  case,  nearly  frequency-independent 
as  we  expect.  I2  is  flat  at  low  frequencies,  and  is  approaching  a  value 


(5.75) 


at  the  left  edge  of  the  graph  in  Fig.  9.  The  limiting  value  of  Ij 
at  the  lower  frequencies,  Ijq  ,  is  related  by  the  ratio  of  acoustic 
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(5.76) 


At  higher  frequencies  becomes  an  increasing  function  of  fre¬ 
quency  as  longer  negligible.  For  example,  at  100  MHz, 

|A22(k2/‘r2)I  =  0.91  for  the  combination  (1  =  methanol,  2  =  aluminum). 

In  the  range  100  MHz  —  1  GHz  the  average  increase  of  Ij  with  w  is 

55  1/2 

«  «'  ,  in  line  with  an  initial  dependence  as  to  ,  approaching 

proportionality  to  cd  at  higher  frequencies. 

Thus  far  we  have  assumed  a  »  T2  ,  and  thus  have  been  dealing 

with  strong  optical  absorptions.  Since  a  is  Independent  of  modulation 

frequency  and  t  «  w  '  ,  at  some  frequency  this  assumption  fails  to  be 

valid  and  v/e  have  rather  different  behavior  than  we  have  seen  thus  far. 

If  we  still  assume  «  Tj,  k2  «  and  «  a  ,  but  we  now 

assume  T2  »  a  ,  and  keep  terms  involving  A2j(k2/a)  ,  then  we  find 

Eq.  (5.66)  and  (5.67)  become: 


-Jk.z  2oI 

a  (z)e  1  - - o_^ pLl 

*  ^21^  L  *"2  J 


1/ 

Ml 

Tj  \ 

1  "2  ' 

\1  +  62i/  « 

(5.77) 


(5.78) 
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Let  us  dispose  of  Eq.  (5.78)  for  a  ,  by  noting  that  since 

•I  "Jk  z 

T2  «  ,  a_j  has  the  form:  a_j(2)e  ^  aiA(l+j)w"^^^  +  B 

where  A  and  B  are  real  constants.  Thus  In  this  limit 

Ij  «  2aV^  +  2AB(ij"^/^  +  . 

The  equation  (5.77)  for  a^g  limit  Is  quite  interesting  in 

2 

that  It  gives  an  acoustic  intensity  proportional  to  a  for  A22  »  1  • 

-5 

Cases  in  point  are  the  combinations  (medium  1  =  10  atm  air;  medium  2  » 

Q 

silicon;  =  6x  10  )  and  (medium  1  =  water;  medium  2  =  silicon; 

^21  '  values  for  A2J  mean  that  both  first  media  look 

more  or  less  like  free  surfaces  to  higher-impedance  silicon.  For  sili- 
con  with  incident  optical  radiation  of  wavelength  8000  A,a=  3.29x10  m  . 
In  this  case  =  1  at  frequency  169  MHz.  For  comparison,  this 

value  of  a  is  a  factor  of  34  smaller  than  that  for  A1  used  in  Fig.  9* 

For  Al,  |T2/»^/a(  *  1  at  483  GHzJ  Hence  no  trace  of  this  effect  will 
be  found  in  Fig.  ‘9,  Equation  (5.66)  is  plotted  as  a  function  of  fre¬ 
quency  for  several  first  media  with  silicon  as  medium  2  in  Fig.  10, — 

In  Fig.  10  we  see  that  for  a  range  of  frequencies  for  the  combina¬ 
tion  water/si 1 Icon  from  about  200  MHz  to  about  2  GHz,  Ig  In  silicon 

Is  a  rapidly  Increasing  function  of  frequency,  the  variation  approaching 
2 

I2  “  w  .  In  this  range  of  frequencies  the  relation  (5.77)  Is  a  good 
approximation  to  the  exact  equation  (5.66),  which  Is  plotted.  For 
example,  for  water/silicon  at  1  GHz,  Eq.  (5.77)  predicts  I-  =  8.7xl0"^® 
W/m  ;  the  exact  value  from  the  graph  1s9.1xio"'®  W/m^. 

The  approximation  In  Eq.  (5.77)  breaks  down  at  the  highest  fre¬ 
quencies  In  Fig.  10  because  k2/o  Is  no  longer  small  and  the  curves 
reach  a  maximum  at  ~  4  GHz. 

-  44  - 


L 


vie  see  that  for  the  approximately  free-surface  case,  in  cases 

where  IT2I  >  a  ,  there  Is  strong  reason  to  prefer  high-frequency 

modulation.  An  immediate  case  in  point  is  electron-beam/acoustic 

microscopy,  where  medium  1  would  be  vacuum,  medium  2  a  solid,  and  a 

would  now  correspond  to  the  electron-beam  attenuation  coefficient.  The 

acoustic  signal  is  detected  with  a  piezoelectric  transducer  in  contact 

2 

with  medium  2.  A  graph  for  vacuum/Si  Ig  would  have  the  same  Ig  *  <0 
characteristic  for  values  of  a  such  that  |T2/v^aI  >  1  .  The  graph 
shown  for  air/silicon  shows  a  signal  increasing  by  two  orders  of  magni¬ 
tude  between  100  MHz  and  1  GHz  for  a  =  3.29  x  lO^m"^  .  A  corresponding 
Increase  will  be  observed  at  lower  frequencies  for  smaller  a. 

To  exhibit  the  dependence  of  acoustic  intensity  on  absorption 
coefficient.  Fig.  11  plots  Eq.  (5.66)  for  Ig  in  silicon  for  the 
combinations  lO"^  atm  air/Si,  water/si,  pyrex/Si,  and  sapphire/Si, 
as  a  function  of  absorption  coefficient  a  .  A  modulation  frequency 
of  100  MHz  is  assumed.  For  high  values  of  a  at  the  right  side  of  the 
graph,  we  find  I2  is  independent  of  a  ,  as  is  predicted  by  Eq.  (5.69). 
In  this  regime  we  get  a  huge  variation  in  generated  acoustic  intensity 
depending  on  the  acoustic  impedance  mismatch  (and  thermal  mismatch) 

across  the  boundary.  Thus  the  I2  signal  level  in  silicon  is  30  dB 

-5 

higher  with  sapphire  as  medium  1  than  with  10  atm  air  (approximating 
vacuum)  as  medium  1.  This  large  dependence  of  acoustic  signal  on  load¬ 
ing  conditions  of  the  surface  has  earlier  been  noted  by  White  and 
von  Gutfeld. 

Physically  we  may  understand  this  in  terms  of  the  reflection  of 
waves  at  the  boundary  between  medium  1  and  medium  2.  If  the  boundary 
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Is  a  free  surface,  the  negative  traveling  acoustic  wave  Is  completely 
reflected  at  the  free  surface  boundary  with  a  180*  phase  change.  When 
a  »  k2  and  iTgl  »  ,  the  positive  traveling  acoustic  wave  and 

reflected  negative  traveling  wave  are  out  of  phase  and  cancel,  leading 
to  small  acoustic  Intensity  1n  the  far  field.  If  medium  1  and  medium  2 
have  equal  acoustic  Impedances,  there  Is  no  reflection  at  the  boundary, 
and  far-field  acoustic  Intensity  comes  only  from  the  positive-traveling 
wave.  If  medium  1  has  infinite  acoustic  impeaance  the  negative  traveling 
acoustic  wave  in  medium  2  Is  reflected  at  the  boundary  with  no  phase 
change.  If  again  a  »  1^2  and  Ixgj  »  kg  ,  the  positive  traveling 
acoustic  wave  and  reflected  wave  are  in  phase,  leading  to  a  large  far- 
field  acoustic  intensity  In  medium  2. 

As  we  move  toward  smaller  values  of  a  In  the  left  side  of 
Fig.  11,  the  graphs  converge  to  near  equality,  as  is  predicted  by 
Eq.  (5.77).  The  water/SI  and  lO'^  atm  alr/Si  curves  In  particular  are 
falling  off  as  in  the  range  of  a  between  lO^m”^  and  10®m'\ 

For  these  materials,  in  the  regime  where  AgjCkg/a)  »  1  ,  Eq.  (5.77) 
simplifies  to 

(5.79) 

or 

'2  “  '» 

2 

which  corresponds  to  the  1/a  dependence  seen  in  Fig.  10  for 
-B 

water/Si  and  10  atm  alr/SI. 


-  47  - 


The  two^medium  casetnaybe  compared  with  White,  whose  classic  paper^ 
considers  absorbtlon  In  an  Infinitely  thin  boundary  layer  of  a  semi¬ 
infinite  medium,  and  their  diffusion  and  expansion  due  to  heating. 

He  considers  two  cases  relevant  to  our  results:  a  perfectly  con¬ 
strained  surface,  and  a  free  surface.  For  the  free  surface.  White 
predicts.  In  the  limit  IT2I  »  k2  *  In  our  notation,  that  far-fleld 
stress 


(5.81) 


where  this  Is  a  rearrangement  of  White’s  Eq.  (5.10).  For  the  case 

vacuum/gold,  Eq.  (5.81)  predicts,  for  f  =  100  MHz  ,  the  far-fleld 

stress  in  Gold  T2  *  13.2x10*®  (N/m^)J  .  Using  Eq.  (5.66)  for  the 

-5  2 

case  STP  air/gold  one  obtains  T2  =  2.6x10  (N/m  )  ,  a  discrepancy 
of  20%  In  stress,  44%  In  power. 

White’s  other  relevant  case,  the  perfectly  constrained  surface, 
predicts  the  stress  In  medium  2  for  the  case  IT2I  »  1^2 

3  V 

IT2I  =  (5.82) 

^2 

where  we  have  rearranged  White’s  Eq.  (5.10)  to  fit  our  notation. 
Recall  that  Eq.  (5.71)  for  far-fleld  stress  In  the  matched  boundary 
case  becomes 


1 

2 


(5.83) 
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when  we  convert  mode  amplitude  a^2  stress.  The  doubling  as  we  go 

from  the  matched  case  (5.83)  to  the  perfectly  constrained  case  (5.82) 

■may  be  attributed  to  the  reflection  of  the  negative-propagating 

acoustics  generated  In  medium  2  at  the  Zqj  =  •  boundary. 

For  example,  In  the  case  sapphire/aluminum  at  100  MHz,  Eq.  (5.66) 

leads  to  stress  Tg  =  1.24xio'^  H/va  whereas  Eq.  (5.82)  gives 
-4  2 

T2  ®  1.73x  10  N/m  .  There  Is  about  3  dB  difference  between  the  two 
values,  somewhat  more  than  in  the  free-surface  predictions  just  con¬ 
sidered.  Considering  the  range  of  acoustic  and  thermal  properties 
available  in  real  materials,  Eq.  (5.82)  seems  somewhat  less  useful  for 
the  case  of  solid-solid  Interfaces  than  Eq.  (5.81)  Is  for  gas-solid  or 
liquid-solid  Interfaces. 

5.  Three-Layer  Model;  Heat  Flow 

We  now  consider  the  3-layer  geometry  of  Fig.  8  In  which  light 
enters  through  transparent  semi -Infinite  medium  1  on  the  left  and  is 
absorbed  In  medium  2,  with  optical  intensity  I(z)  =  In  medium  2. 

Medium  2  Is  an  Infinite  slab  of  thickness  i  which  contacts  medium  1 
at  2  -  0  and  semi-infinite  medium  3  at  z  =  t  .  Medium  3  Is  assumed 
to  be  transparent.  Our  procedure  Is  to  first  solve  for  the  temperature 
.butlon  In  space,  and  then  use  this  as  the  distributed  source  for 
positive  and  negative  traveling  acoustic  modes.  We  shall  be  concerned 
finally  with  the  far-field  modes  in  media  1  and  3,  and  their  associated 
acoustic  power.  Parker^®  studied  the  case  where  medium  1  was  glass, 
medium  2  a  negligible-thickness  absorbing  layer,  and  3  the  gas  In  a 
photoacoustic  cell.  Rosencwaig  and  Gersho  solved  for  the  pressure  In  a 
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I 


gaseous  medium  1  when  medium  2  was  an  absorbing  solid  and  medium  3,  a 

22 

backing  material,  von  Gutfeld  and  Budd  considered  a  3-layer  case  where 
medium  2  absorbed  all  optical  power  In  an  Infinitesimally  thin  surface 
layer.  They  found  the  acoustic  signal  In  liquid  medium  3  to  be  strongly 
dependent  on  layer  thickness  and  choice  of  liquid. 

Our  analysis  to  follow  is  applicable  to  any  choice  of  media  1,  2, 

3  and  we  arrive  at  exact  analytic  expressions  for  far-fleld  acoustic 
amplitudes  with  no  Initial  assumptions  on  the  relative  sizes  of  optical 
and  thermal  skin  depths  or  acoustic  k-vector. 

We  proceed  as  In  the  single-interface  case.  In  medium  2,  pro¬ 
ceeding  as  In  Eqs.  (5.16)  and  (5.17),  but  noting  that  Integration  In 
the  latter  Is  now  from  z  to  i  ,  we  arrive  at  the  3-layer  analogs  of 
Eqs.  (5.19)  and  (5,20): 


(<• 


al  r  -T«Z  -azl  X  -T^z 

Ttil' 


(5.84) 


=  -2 - - 

(a  +  T«)  L 


a_2(Jt)e 


(5.85) 


Since  the  heat  source  exists  only  in  medium  2,  In  media  1  and  3  we  will 
have: 


a^j(z) 


a^,3(z) 


3:^(0)  e  1 

(5.86) 

.  -Tx(z-i) 

4(1)6  3 

(5.87) 
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and  aJjCz)  =  ~  ^  boundary  conditions  are  that  heat 

flux  f  and  temperature  6  are  continuous  at  z  =  0  and  at  z  ~  I  . 
Application  of  these  permits  us  to  evaluate  the  unknowns  a^2(0}, 
42(0),  .^2  (t)  and  •  Then  in  region  i  we  recall  that 

temperature  0^(z)  =  (a^^(z)  -  a^. (z))/22^  .  Solving  for  9,.(z)  in 
region  1,  we  have: 

T,Z 

0l(z)  =  Gj^e  ^  (5.88) 


In  region  3: 


03(1)  -  6356 


~T,(Z-Jl) 


(5.89) 


In  region  2: 


-Tj,Z  +Tj,Z  -az 

62(2)  -  ®20^  ^ 


'21' 


'22' 


(5.90) 


where  solution  of  the  boundary  condition  equations  yields: 


«  _  ^21®02 

®11 - 

Dt 


T2A 


®23®02 


'30 


-?r) 


(a-Tg)^ 


(5.92) 


*20 


Dt  L 


-  (1  -  623)  (1  -  ^23  *2  )®  j 


(5.93) 


»21  “  ^  ■  '23^  (1  "  «21  ^  «2l)  (1  -  «23  ^) 


-al 


(5.94) 


«22  '  -®02 


(5.95) 


where 


«tj  '  A'^l 


eQ2  given  by  Eq.  (5.31),  and 


(5.96) 


Dt  *=  (1  -  621)^1  ■  <*23^®  ^  ^23^®^^^ 


(5.97) 


These  solutions  for  the  temperature  distribution  in  the  3-medium 

geometry  agree  with  the  harmonic  part  of  the  solution  given  by 

21 

Rosencwaig  and  Gersho  for  a  similar  3-medium  geometry. 

6.  Three-Layer  Model:  Acoustic  Generation 

With  the  temperature  distribution  in  the  3  regions  known,  we  can 
solve  for  the  acoustic  modes.  Since  there  are  sources  [non-vanishing 
e^(z)l  in  all  3  regions,  we  will  have  positive  and  negative  traveling 
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inodes  In  each  region.  For  region  1, 

. 


Z 

(z)  =  ^  ^  y*  d;  +  a^j(-«)e  ^  (5.98) 


where  the  third  subscript  on  c^jj  designates  the  region.  Since 
*+1^““^  =  0  (no  acoustic  sources  at  “®  )  with  Eq.  (5.88)  we  are 
led  to: 


JkjCjjjSjOjj  TjZ 


Tj  t  jk, 


(5.99) 


Similarly,  for  the  negative  traveling  mode 


jk.z  -0  -jk.i;  jk.z 

.l(z)  =  ^  /  0j(;)e  ^  dt  +  a_j(0)e  ^ 

•4 


(5.100) 


which  yields 


jk.c..,8.e,,  jk.z 
a  ,(2)  =  -IJll.l  .11  e  1 

(jkj-Tj) 


-(jki“T,)z 

1  -  e  ^  ^ 


jk.z 

+  a_j(0)e  ^  (5.101) 


In  region  3  we  have: 


jk-z  /•  -jkt 

a_3(z)  =  “jk3Cjj3P3e  ^  j  03(c)e  H  +  a_3(-H»)e 

•z 


jk-z 

^  (5.102) 


No  acoustic  sources  at  z  =  +«  means  a_3(^)  *  0  ,  and  then  substi¬ 
tution  of  Eq.  (5.89)  results  In: 


“jk^jC.  .  ”''^o(z*'J') 

a  ,(z)  =  — IMIUO  e  3  ^ 

-3  (jkj  *  Tj) 


-T.(Z-JI) 

3  (5.103) 
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For  the  positive- traveling  mode  in  region  3  we  will  have: 


•jk  2  *  jk  t 


-jkjCz-il) 


(5.104) 


We  substitute  Eq.  (5.89)  for  0^  and  get: 


a+3(2) 


~'^S^113^3^30 

(jka  -  T3) 


-T-(2-t)  -jk,(z-t) 

e  -  e 


+  a^3(t)e 


-jk3(z-t) 


(5.105) 


In  region  2  we  will  have: 


,  ~jk-z  f  jkpt  -jk^z 

*+2(2)  “  -JV112V  /  +  i+2^0)e  ^  (5.106)- 


Substitution  of  Eq.  (5.90)  in  this  expression  yields; 


-jk-z 

a^.2(2)  »  -jk2Cjj202® 


®20_  /  _  j\  ^  ”21 

To  '  '  jk2  +  T2 


l_  jk2  ‘2 


/  (jkp+T„)z  V  Opp  /  (jkp-a)z  v' 

r  -f 


(0)e 


-jk^z 


(5.107) 

Let  us  define  E22  as  the  entire  first  term  on  the  right  evaluated 
at  z  =  t  ,  that  is: 


®+2^^^  “  ^22  ® 


■jk„t 


(5.108) 
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Likewise,  the  negative  mode  will  be: 


= 


jkgZ 


X. 

/ 


,({)e  ^ 


d5  + 


a_2U)e 


jk,(z-A) 


(5.109) 


Substitution  of  Eq.  (5.90)  and  Integration  gives  us: 


2  ■  ‘2^ 

022  /  “(ok2“T2)t  ”(jk^“T(,)2 

,-tJ  V 


(jk2-T2) 


2“i2^*  \j''2~^2^^^ 


®22  / 

(jk2  +  a)  V  ' 


+  a^2^^)® 


jk2(z-t) 


(5.110) 


We  define  E2Q  as  the  entire  first  term  on  the  right  side  of  this 
equation  at  z  =  0  ;  that  is: 

-jk  1 

a_2(0)  =  E20  +  a_2{t)e  ^  (5.111) 

The  remaining  4  boundary  conditions  (the  first  2  were  vanishing 
a^j(-“)  and  vanishing  a_2(+»))  are  equality  of  stress  T  =  (a^+a_)/2 
and  particle  velocity  u  =  (a^-a_)/2Zp  on  both  sides  of  the  boundaries 
at  z  «  0  and  z  =>  i  . 


i 

i 
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T  at  z  »  0  : 

a+j(0)  +  a^j(O)  « 

•  jk  fi 

or  Ej  +  a^j(O)  -  a^2t®)  ^20  ^ 

a^gCO)  -  a_2(0) 

-jkgi 

a^2^®)  “  a_2(t)6  “  EgQ 

(5.113) 

T  at  z  *  1  : 

-jkoJl 

or:  *+3^*-)  ■’■  ^3  *  a^2^®)®  ^22  *  ®-2^*'^  (5.114) 

V  at  z  *=  A  : 

^23[®+3^^^  ‘  ^3]  '  ®+2^°^®  ^  ^22  ■  ®-2^*'^  (5.115) 

Equations  (5.112),  (5.113),  (5.114).  and  (5.115)  are  a  set  of 
4  equations  in  the  remaining  4  unknowns:  a  j(0),  a^g^^)*  a_2(t)» 
and  a+3(A)  .  We  need  find  only  a^j(O)  and  8^3(4)  to  write  down 


u  at  z  «  0  : 


|52  -  >_j(0)]  . 

^01 

or:  ^2l["®-l^®^  ^  ^ll  * 
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the  far  field  behavior.  Solution  yields: 


^  [(‘  ^  ^  +  (1  +  A23)(1  -  A2j)e  ^  jE^ 


jk  1 

-  2(1+A23)e  ^  E2q“2(1-A23)E22-4A23E2 


»+3(t)  =  -  jj  *  ^^^"^21^^20  *  ®  ^ 

‘•Jk  A 

+  ((1-A2i)(l+A23)e  ^  -  (1  +  ^2^(1  - «23> 


(5.116) 


jkpi.  I 

-  (1  +A2j)(1  - A23)e  jE3j 

(6.117) 


where 


“jk  t  jk  t 

0  *=  (l“A2i)(l-A23)e  ^  -  (1  +  A2j)(l  +  A23)e  ^  (5.118) 


Letting  z  approach  “*  in  Eq.  (5.101)  and  letting  z  approach 
in  Eq.  (5.105)  results  in  the  far  field  modes: 


jk. c. . . 8.6. .  jk.z  J 

a  (z)  »  -1-1.11..U1  e  i  +a,(0)e 

(jkj-Tj) 

a  (z)  =  ^*^3*^113^3^30 


(5.119) 


(jk3  -  T3) 


-jk,(z-i) 

+  a^3(t)e 


(5.120) 


If  we  now  substitute  (5.116)  and  (5.117)  in  these  equations,  making 


-  57  - 


use  of  (5.103)  and  (5.99),  we  obtain; 


2  ( 


-  (1  +  A23)(jk 


+  (1  +  A23)e  E2Q  +  "A23)E22'''^23^3 


jkpAl 

Ci+A21^iJ®  J 


*«<*)  •  -5 


E,  r  -jk.t 

(jka-ta)  ‘^21^^'^*^3‘^23''3^®  -(1+A2j) 


\5.121) 


>«  (jk^  +  A23T3)e 


jk2J. 


jkp"  )  -jk,(z-t) 

+  (l-A2i)E2Q+(l+A2j)e  E22  +  ^AgjEj  }  e 


(5.122) 


These  expressions  are  our  final  equations  for  the  far-field  acoustic 
inodes  in  regions  1  and  3.  Since  a^j(2)  =0  in  the  far  field  of  region  1, 
and  8^3(2)  ®  0  in  the  far  field  of  region  3,  these  expressions  give  the 
far-field  stress  and  particle  velocity,  and  acoustic  power: 


(far-field  stress) 

(far-field  particle  velocity) 


1^(2)  » 


a_j(2) 


Uj(2) 


a_l(2) 


2Z 


01 


(5.123) 

(5.124) 
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(far-field  acoustic  intensity) 

'^01 

(5.125) 

(far-field  stress) 

^  2 

(5.126) 

(far-field  particle  velocity) 

.3(z)  = 

^  21 
^^03 

(5.127) 

(far-field  acoustic  intensity) 

1  a^3(z)a  -(z) 

I  (z)  =  - +3 — 

8 

(5.128) 

We  now  consider  the  examples  listed  in  Table  5.3.  A  Fortran 
program  to  calculate  far-field  acoustic  intensities  and 
from  Eqs.  (5.125)  and  (5.128),  using  Eqs.  (5.121)  and  (5.122)  was 
implemented.  This  program's  listing  is  included  as  Appendix  D. 

First  we  consider  the  effects  of  the  third  medium,  which  is  the 
new  element.  The  effects  of  adding  the  third  medium  on  both  Ij  and 
can  be  profound.  If  medium  3  is  an  efficient  thermal -acoustic 
generator,  such  as  acetone,  large  signal  enhancements  are  possible, 
for  very  thin  films,  as  von  Gutfeld  and  Budd  have  observed.  Also, 
the  2-3  boundary  means  we  can  get  both  a  reflected  wave  in  medium  2, 
and  thickness  resonance  effects. 

In  Fig.  12  and  Fig.  13  we  see  intensity  1^  and  I^  in 

the  far  field  of  medium  3  for  pyrex  as  medium  1,  aluminum  as  medium  2, 
and  4  different  third  media:  acetone,  water,  methanol,  and  aluminum. 
Aluminum  as  the  third  medium  reduces  the  3-layer  problem  to  the  2-layer 
case  already  studied  and  is  included  for  comparison. 
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TABLE  5.3.  Photoacoustic  Generation  In  3-reg1on  Systems 
Discussed  In  Section  5.7 


Graphs 


Figure  Points  Illustrated 


Pyrex/Al /Acetone 
Pyrex/Al/Al 

Pyrex/Al /Methanol 
Pyrex/Al /Water 

5-12 

Effect  on  Ij  and  I 3  of  changing 
medium  3  for  different  layer 
thickness. 

Pyrex/Mo/ Acetone 

Sapphi re/Mo/Acetone 

5-14 

Effect  on  of  changing 

medium  1  for  different  layer 
thickness. 

Pyrex/Al /Water 

Py r ex/ Au/ Water 

Pyrex/SI /Water 
Pyrex/Mo/Water 

5-15 

Effect  on  I3  of  changing 
medium  2  for  different  layer 
thickness. 

Pyrex/Mo/Water 

5-16 

Effect  on  I3  of  changing 
modulation  frequency  for 
different  layer  thickness. 

Sapphi re/Carbon/Acetone 

5-17 

Effect  on  I,  and  I,  of  changing 

Sapphi re/Carbon/Water 

5-18 

layer  thickness.  A/4  thickness 

resonance. 

Sapphlre/Al /Water 
Sapphire/Si /Water 

5-19 

Contrast  in  a  photoacoustic 
image. 

W/m^  I|.  W/m 


PYREX/Al  WITH  VARIOUS  MEDIA  3 
I|  IN  PYREX  lOpMHz 


L- 1.0/1 

NORMALIZED  LAYER  THICKNESS  A2L 
FIGURE  12 


PYREX/Afi  WITH  VARIOUS  MEDIA  3 


L«  I.O/i 


NOmALiZED  LAYER  THICKNESS  Ag  L 
FIGURE  13 
-  61  . 


For  the  case  of  aluminum  at  an  optical  wavelength  0.95y  ,  optical 

Q 

absorption  constant  Is  a  =  1.12  x  lo  m  .  For  100  MHz  modulation, 
thermal  attenuation  constant  82  *  1^2/*^!  =  1*61  x  10®  m‘^  .  The 
two  graphs  are  plotted  as  a  function  of  normalized  layer  thick¬ 
ness  82!-  .  The  effects  of  finite  {rather  than  Infinite)  a  are  seen 
only  at  the  extreme  left,  where  aL  s  1  or  less,  and  thus  a2L  «  1  . 

In  this  region  all  the  curves  approach  zero  as  a2L  -*•  0  ,  since  not 
all  the  light  Is  being  absorbed. 

For  the  case  of  In  the  region  from  roughly  a2L  =0.03  to 
82!-  ®  0.3  (200Ato  2000A  thickness)  In  Fig.  13.  both  methanol  and 
acetone  give  a  large  signal  enhancement  over  the  2-reg1on  case  where 
the  A1  extends  to  +•  .  The  Intrinsic  conversion  factors  n  from 
Table  5.2  are  In  the  ratio 

Acetone :Methanol :A1 :Water;:  185:81 :2;1 

The  Intensities  In  Fig.  13.  are  In  the  same  order  for  thin  A1  films. 
Physically,  since  a2L  «  1  ,  most  of  the  heat  generated  at  the  pyrex/ 
aluminum  Interface  reaches  the  aluminum/medium  3  Interface,  where  the 
transmitted  heat  wave  In  the  third  medium  generates  acoustics  In 
medium  3  by  the  expansion  of  medium  3.  When  medium  3  has  a  large 
Intrinsic  conversion  factor  n  ,  large  signal  enhancement  Is  possible. 

For  thicker  layers  we  have  a  diminishing  amplitude  thermal  wave 
reaching  medium  3,  and  the  acoustic  wave  Intensity  decreases  as  the 
square  of  the  thermal  wave  amplitude.  For  a2L  >2  In  the  present 
case,  essentially  all  the  sound  reaching  the  far  field  of  region  3  Is 
created  In  region  2  (Al)  and  transmitted  across  the  boundary  from  2  to  3, 
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When  medium  3  is  one  of  the  liquids  the  impedance  mismatch  at  this  bound¬ 
ary  with  A1  is  large.  For  e^cample,  ^23  =  =  18.7  for  medium  2  = 

Al,  medium  3  *  acetone. 

For  this  reason,  the  signals  I3  in  the  liquids  are  considerably 
smaller  for  thermally  thick  films  than  for  the  matched  case  when 
medium  3  is  Al.  We  note  that  if  one  were  to  insert  a  fourth  medium 

^04  ~  *^02^03  thickness  between  medium  2  and  3,  the 

signals  in  the  liquids  would  rise  to  the  same  value  as  the  2-medium 
pyrex/Al  case,  as  shown  in  Fig.  12. 

Resonance  behavior  is  not  much  in  evidence  for  the  I3  graphs 
only  because  medium  1,  in  this  case  pyrex,  is  a  rather  good  acoustic 
match  to  Aluminum  ~  '^21  *  1*33)  and  thus  the  power  reflection 

coefficient  of  a  wave  traveling  toward  medium  1  in  medium  2  is  only 
0.061.  We  will  see  examples  later  where  I3  exhibits  strong  resonance 
effects  with  substantial  mismatch  in  impedance  between  media  1  and  2. 

Looking  at  the  far-field  acoustic  intensity  in  pyrex  for  the 
same  four  combinations,  Fig.  12,  we  find  that  in  every  case  the 
addition  of  a  third  medium  decreases  the  signal  for  the  range  of 
thickness  shown.  The  sharp  increase  in  as  aL  -»■  0  for  the 
methanol  and  acetone  graphs  is  to  be  understood  as  sound  produced  in 
region  3.  The  liquids  with  high  conversion  factors  n  receive  a  heat 

wave  produced  in  medium  2  (Al)  near  the  1-2  boundary  which  has  not 
been  completely  damped  by  the  time  it  reaches  the  2-3  boundary,  for 
a^L  ;$  2  .  Sound  produced  in  medium  3  propagates  in  the  negative  z 
direction  as  well  as  the  positive  z  direction.  A  percentage  (19%  in 
the  case  of  medium  2  =  aluminum,  medium  3  =  acetone)  of  the  negative¬ 
traveling  sound  power  in  medium  3  crosses  the  2-3  boundary  into 
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medium  2,  where,  because  of  the  good  Impedance  match  between 
medium  1  (pyrex)  and  medium  2  (aluminum),  there  is  negligible  reflec¬ 
tion  at  the  1-2  boundary,  and  the  sound  travels  to  the  far  field  of 
medium  1.  This  accounts  for  the  sharp  rise  in  Ij  for  a  very  thin 
region  2. 

The  Ij  graphs  in  Fig.  12  show  a  drop  of  almost  20  dB  from  the 
2-medium  case  of  pyrex/Al  for  thermally  thick  (but  still  acoustically 
thin)  layers  near  aL  ~  1.5  .  To  explain  this,  we  should  note  that 
in  no  case  is  much  sound  produced  in  pyrex  itself.  From  Table  5.2, 
^2^’^1  *  ^  2  =  aluminum.  The  substantial  signal 

in  the  2-layer  pyrex/Al  case  is  sound  generated  by  thermal  expansion 
of  A1  which  crosses  the  1-2  boundary.  When  we  have  a  third  low- 
impedance  medium  terminating  region  2  for  L  «  X  ,  the  wave  traveling 
toward  +2  in  medium  2  is  reflected  at  nearly  full  amplitude  but  at 
opposite  phase.  This  wave  Interferes  with  the  wave  traveling  towards 
“2  originally,  resulting  in  drastically  reduced  Ij  .  As  the  layer 
becomes  thicker,  increases  as  the  two  waves  get  back  in  phase 
with  each  other,  which  happens  at  the  X/4  thickness,  a^L  s  25  for 
100  MHz,  not  shown  in  Fig,  12.  This  would  be  a  16  p  thick  A1  region. 

To  see  the  effect  of  the  first  medium  on  the  sound  transmitted 
to  the  far  field  of  medium  3,  we  compare  the  cases  Sapphire/Mo/Acetone 
with  Pyrex/Mo/Acetone  in  Fig.  14.  The  for  Sapphire  is  lower 
for  thin  film  thickness  because  Sapphire  is  a  better  heat  conductor 

than  pyrex  ^^sapph1re'*^*^pyrex  *  reach 

acetone,  whose  very  efficient  expansion  is  responsible  for  the  far- 
fleld  sound  in  the  case  of  a  thin  Mo  film.  In  the  thick-film  case,  the 
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SAPPHIRE/MO//ACETONE  I3  100 MHz 
PYREX/M0//ACET0NE  I,  100  MHz 


NORMALIZED  THICKNESS  A 

FIGURE  14 


sound  reaching  the  far  field  of  medium  3  is  actually  produced  in  media  1 

and  2.  Sapphire  has  a  much  higher  intrinsic  conversion  factor  n  in 

Table  5.2:  ®  710  .  It  is  therefore  more  efficient  as 

sapphire  pyrex 

a  photoacoustic  generator  than  pyrex,  and  is  in  fact  more  efficient  than 
molybdenum.  It  is  also  a  more  efficient  constraint  on  molybdenum 

^^o.sapphire  ®  "'o.pyrex^* 

field  acoustic  power  in  acetone  in  the  thick  film  case. 

We  next  consider  the  effects  of  a  change  in  the  second  medium  when 
the  first  and  third  media  are  kept  fixed.  In  -Fig.  15  we  have  again 
plotted  Eqs.  (5.128),  using  Eq.  (5.122)  for  four  different  second  media: 
gold,  molybdenum,  aluminum,  and  silicon.  Medium  1  is  pyrex  in  all  cases 
and  medium  3  is  water.  There  are  several  factors  to  consider  to  account 
for  the  differences  we  observe. 

Firstly,  the  silicon  curve  has  a  more  gradual  rise  at  small  film 
thickness.  This  is  because  the  optical  absorbtion  constant  a  for 
silicon  (3.29x10®  m"^,  optical  X  =  8000  A)  is  a  factor  of  at  least  13 
smaller  than  the  a  values  for  the  other  materials.  A  thicker  Si  layer 
is  required  to  absorb  all  the  light. 

The  curves  for  A1  and  Si  in  Fig,  15  are  flat  as  a  function  of 
layer  thickness,  whereas  the  Mo  and  Au  graphs  are  decreasing  functions. 
The  difference  is  the  acoustic  impedance  match  with  medium  1,  pyrex. 

Aluminum  (A22  *  ^02^^01  '  (^21  "  fairly  well 

matched  to  pyrex.  For  a  wave  in  medium  2  striking  the  1-2  boundary, 
the  power  reflection  coefficients  are  0.02  and  0.04,  respectively,  versus 
0,47  and  0.42  for  Mo  and  Au.  This  means  that  acoustic  thickness  reso¬ 
nance  effects  will  not  be  much  in  evidence  in  the  far-field  in  medium  3, 
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(  IN  WATER  FOR  VARIOUS  MEDIA  #2 
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there  being  only  the  portion  of  the  positive  traveling  wave  In  medium  2 
that  is  transmitted  across  the  2-3  boundary  to  consider.  Gold  and  Mo, 
being  poor  acoustic  matches  to  pyrex,  exhibit  thickness  resonance  effects 
due  to  multiple  reflections  within  the  layer.  Both  exhibit  a  half-wave 
resonance  maximum  (at  thickness  34pfor  Mo,  16pfor  Au)  and  a  quarter- 

wave  minimum.  In  the  thin  layer  region  shown,  the  Mo  and  Au  curves  are 
decreasing  towards  their  quarter-wave  minimum  off  the  right  side  of 

the  graph.  The  details  of  the  Mo  layer  thickness  resonance  will  be 
considered  shortly.  We  emphasize  that  because  of  the  large  acoustic 
Impedance  mismatch  between  all  four  solids  and  water,  the  sound 
intensity  In  water  Is  very  low  compared  to  the  sound  Intensity  we 
would  have  If  medium  2  had  extended  to  “  . . 

In  Fig.  16  we  have  plotted  Intensity  In  the  far  field  In  water 
for  the  combination  pyrex/molybdenum/water,  also  studied  by  von  Gutfeld 
and  Budd.  We  see  here  the  effect  of  changing  modulation  frequency  of 
the  Incoming  light.  In  the  region  at  the  extreme  left,  for  film  thick¬ 
ness  approaching  zero,  the  far-fleld  acoustic  Intensity  goes  to  zero 
because  the  film  thickness  Is  less  than  the  optical  penetration  depth 
and  the  light  Is  not  all  absorbed.  For  optical  X  =  1.0  y  ,  In 
molybdenum  1/a  =  224  A  ,  In  accord  with  the  turning  point  at  the 
extreme  left  of  each  graph  in  Fig.  16. 

In  the  region  of  film  thickness  from  roughtly  0.04  y  to  ly  ,  the 
far-fleld  Intensity  I^  In  water  decreases  with  increasing  film  thick¬ 
ness  at  each  frequency.  This  Is  because  the  thermal  wave  in  the  molyb¬ 
denum  film  Is  attenuated  with  attenuation  constant 
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FIGURE  16 


Thus  a  «  /jj  means  a  smaller  amplitude  thermal  wave  reaches  the  water. 
Water  has  approximately  the  same  Intrinsic  conversion  factor  n  as 
molybdenum,  as  we  see  In  Table  5.2.  Water  Is  a  better  thermal  match 
to  molybdenum  q  -  9*6)  than  an  acoustic  match 

Z..  u  n  “  45)  .  Thus  sound  In  the  far  field  of  water  Is  most  effi- 
ciently  generated  by  the  thermal  wave  heating  the  water  directly, 
since  sound  generated  Inside  the  Mo  film  doesn't  get  out  efficiently. 

This  analysis  breaks  down  In  the  right  hand  region  of  the  1  GHz 
graph  In  Fig.  16.  At  film  thickness  3.4  u  the  1  GHz  graph  reaches 
a  peak  about  3  dB  below  the  very  thin  film  peak  at  about  0.04  y  thick¬ 
ness.  At  3.4  y  thickness  the  film  acts  as  a  half-wave  acoustic 
resonator,  since  the  acoustic  wavelength  in  Mo  at  1  GHz  is  6.8  y. 

Due  to  the  large  reflection  coefficients  at  each  boundary,  acoustic 
waves  in  the  Mo  film  are  reflected  many  times  back  and  forth  within 
the  layer,  some  energy  being  transmitted  into  water  on  each  pass. 

The  transmitted  waves  in  water  are  all  in  phase  when  L  =  \/Z  ,  and 
hence  we  get  maximum  intensity  with  this  condition,  as  Fig.  16  shows. 
The  A/2  thickness  resonance  condition  is  appropriate  when  both  bound¬ 
aries  are  low  impedance  boundaries  as  in  this  case. 

Cases  exhibiting  a  quarter-wave  resonance  thickness  for  the  cen¬ 
tral  layer  are  sapphire/carbon/acetone  and  sapphire/carbon/water,  for 
»di1ch  and  Ij  are  graphed  as  a  function  of  film  thickness  in 
Figs.  17  and  18.  The  carbon  data  are  an  average  of  thin-film 
carbon  specimens.  The  thermal  diffusion  constant  a  for  carbon  at  . 

1  GHz  is  such  that  aL  =  7.5  corresponds  to  film  thickness  L  =  l.Oy  . 
The  sound  wavelength  in  carbon  at  1  GHz  is  4y  ,  giving  us  the  quarter- 
wave  resonator  we  see. 


W/m^  I3,  W/m 


We  expect  quarter-wave  resonance  In  this  case  because  carbon 
Zq  ■  8x10®)  Is  intermediate  in  acoustic  impedance  between  sapphire 
(z^  “  44x10®)  and  water  (z^  =  1.5x6)  .  In  general,  L  =  , 

where  is  an  o^  integer,  satisfies  this  resonance  condition.  In 
fact,  carbon  is  very  close  to  the  geometric  mean  (z  =  8.2 x  lo®)  between 
the  impedances  of  sapphire  and  water  and  thus  acts  as  an  ideal  quarter- 
wave  transformer  between  sapphire.  Thus  the  acoustic  mode  traveling  to 
the  right  in  sapphire  (medium  1)  in  Fig.  8  is  not  reflected  at  the  ' 
sapphire/Mo  boundary.  Similarly  the  acoustic  mode  traveling  toward 
medium  2  in  medium  3  is  not  reflected.  In  this  matching-layer-resonator 
case  we  expect  Ij  and  I^  to  be  equal.  We  see  from  Figs.  17  and 
18  that  Ij  and  are  very  nearly  equal  for  L  =  X/4  s  7.5/32 
for  Sapphire/Carbon/Water.  This  is  somewhat  less  true  for  Sapphire/ 
Carbon/Acetone,  where  the  ideal  matching  layer  would  have  z^  *  6.4x10® 
mks  rayl,  20%  less  than  the  carbon  value. 

The  resonance  behavior  of  photoacoustic  generation  in  layers  is 
important  in  that  layer  thickness  is  then  an  important  source  of  con¬ 
trast  in  photoacoustic  images.  One  can  imaging  extracting  layer  thick¬ 
ness  information  from  photoacoustic  Images  if  one  can  tune  modulation 
frequency  over  a  substantial  range. 

Finally,  to  illustrate  the  application  of  these  results  to  a  real 
photoacoustic  image,  we  have  plotted  the  response  I^  in  water  for  the 
cases  Sapphire/Aluminum/Water  and  Sapphire/Silicon/Water  in  Fig.  19. 

In  Figs.  5  and  6  we  saw  areas  of  this  nature  and  noted  the  silicon  on 
sapphire  areas  gave  much  weaker  response  than  the  aluminum  on  sapphire 
areas.  The  aluminum  areas  are  1.0  p  thick  while  the  polycrystalline 
silicon  regions  are  0.5  p  thick.  Figure  19  predicts  a  difference  in 
response  in  the  far  field  of  ~  11  dB,  consistent  with  what  we  saw  in 
Figs.  5  and  6. 


SAPPHIRE /Aje/WATER  I3  840MH2 
SAPPHIRE /Si/ WATER  I,  840  MHz 
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VI.  OTHER  EVENTS 


As  mentioned  earlier  we  have  had  visits  from  Dr,  Milton  Ritchie  and 
Dr.  Oim  Knight  from  the  Hughes  Santa  Barbara  Research  Center.  They  are 
asking  us  to  study  the  adhesion  of  lead  sulphide  films  in  the  infrared 
detectors  —  a  device  that  is  being  built  under  contract  with  the  Air  Force 
(Space  Division). 

Dr.  J.  Reffner  of  American  Cyanamide  Company  has  supplied  us  with  a 
number  of  samples  of  polymers  for  study  with  acoustic  microscopy. 

A  group  from  Xerox  Research  Laboratory  in  Toronto  has  visited  us  and 
they  left  with  a  strong  interest  in  our  work. 

Jack  Baring  of  Honeywell  has  supplied  us  with  samples  of  magnetic 
recording  heads  that  are  defective  in  the  sense  that  they  generate  exces¬ 
sive  second  harmonic. 


List  of  Publications 

“Scanned  Acoustic  Microscopy,"  E.  A,  Ash,  ed..  Chapter  II,  Academic  Press, 
1980. 

"Film  Adhesion  Studies  with  the  Acoustic  Microscope,"  by  R.  C,  Bray, 

C.  F.  Quate,  0.  Calhoun,  and  R.  Koch;  published  in  Thin  Solid  Films, 
74,  295-302  (1980). 

Seminars 

“Microwaves  and  Acoustic  Imaging,"  Physics  Department  Colloquim,  Princeton 
University,  Princeton,  PA 

"Acoustic  Imaging,"  Electrical  Engineering  Department,  Stanford  University. 
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Some  of  our  Visitors 


T.  Oestephano  and  V.  Jipson  of  IBM  at  Yorktown  Heights,  New  York 

Or.  Akiyama  of  Hitachi  in  Japan 

Dr.  J.  Zelich  of  Hughes  Aircraft  Company. 

Dr.  P.  Liao  of  Bell  Labs 

Dr.  H.  Fritsche  of  Leitz  in  West  Germany 

Dr.  K.  Wishnuff  of  National  Semiconductor 

Dr.  D.  Jonney  of  Los  Alamos,  New  Mexico 

Dr.  J.  McGroddy  of  IBM  at  Yorktown  Heights,  New  York 

Dr.  C.  Tracy  of  GM  Research  Labs. 

Dr.  F.  Jamerson  of  GM  Research  Labs. 

Dr.  K.  McKay  of  Avantek. 
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PBOCRAH  TO  CklCOLAlE  THE  BEP1ECT4SCE  EDDCTIOB 
OP  H  LAXEFS  CM  «  SUEST8ATE  AS  A  PCHCTIOH  OP 
SIVE  OP  ANGIE  Cl  IKCIOEHCE. 


//BEPLEC  JOB  'f  17XE1.1S1*«<B  BBAt'«CLASS>I 
/♦JOEPARn  DES1=SIIE 
//  EXEC  WATPIV 
//GO.SYSIH  DO  • 

SUATPIV 

C 

C  THIS  PBOGBAE  CAICULATES  TEE  BEPLECTAKCE  FDKCTION  (HAGBITOCE 
C  ABO  PHASE)  FOB  A  ElAAE  ACOUSTIC  NIVE  INCICEMT  ON  A  RULTI- 
C  LAYER  TARGET.  AS  A  FURTION  CP  AAGIE  CF  IHCICENCE.  IT  ALSO  CALCULATES. 

C  TRAMSHISSICN  FECK  THE  LICUIC  TEFOOCE  KANY  lAIERS  INTO  A  SOLID  . 

C  THE  HATHEMATICS  PCILCWS  RA  lEHCAS*  THESIS,  "ACOUSTIC  HICROSCOPY 
C  BY  MECHANICAL  SCANNING",  STASFCSO  ONIVERSIII,  J975,  BUT  DIFPEBS 
C  IN  THAT  THE  INCIDENT  HAVE  IS  IN  THE  LIQUID.  ISOTROPIC  SOLIDS 
C  ARE  ASSUMED.  . 

C 

COMflOH/CCAIC/AlPHA (8) , BETA (8) ,V  |B) , FP (8) , C (8) ,0P(8), 

1C11  (8)  ,C12  (8)  ,C4(J  (8)  ,K  <8)  .KAPPA  (8)  , SIGMA, NP2 
REAL*8  THICK  )9),EBO(9) ,C  11, C 12, C4«, SIGMA, K, KAPPA, 

2HAGAOP(4  6,2)  ,  M AGECP  (4  6 ,2) , MAG  A3  (4  6,2)  .THETA, SIGSQ,TH, PI, Xl,X2,X3, 
3X4,X5,X6, EHTACF (46,2) , P RIECF (46 ,2) ,PH1A3(46,2) ,COSTU, 
4XXR,YXI,YXr,XXF,TLM(100) ,TLF (ICC) ,TEM  (100) , ISP  (100) , 

5RLK(100)  ,KIF  (1C0) 

COBPLEX*10  Z (3,3) ,V(3) ,t (3,3) ,DFT, DET 1 ,DET2, DET3, AOP, BOP, A3,DCHPLX 

1, ANPl.BNF1,A0ER,JAr,PRCD(4,4) , AIPHA.EETA 

2, nATl (4,4)  ,  MAT  2  (4,4) ,ICIRA1(4,4} ,F,PP,Q,QP 

PI"3.1415926S35eS79 
JAY=(0.0CC,1.0C0) 

C 

SSPACE  3 

READ,H 

NP1=N+1  -  . 

NP2"Nt2 
KH1=N-1 
DO  701  1*1, B 
C11  (I)*1.0 
C44(I)  *1.C 
701  RHO(I)  =  1.0 

DO  100  1*1, NP2 

100  R£A0,C11(I),C44(I),BN0(1),THICK(I) 

DO  no  1=1, HP2 
C12(I)  =C11  (I)-2.»C44  (T) 

110  K(I)=0SCFT<PHC(I)/C11(I)) 

DO  115  1=2, KP2 

115  KAPPA(I)  =  CSCR1  (BH0(I)/C44(I)) 

HRITE  (6,50) 

50  P0RMAT(1II1,  12X,»S01I0  C',5X,*SCLID  1',5X, 'SOLID  2',SX,'SCI1D  3', 
ISX.'SOllC  «• ,5X, 'llOUIt',/) 

HRITE  (6,51)  (C 11 (I) ,1  =  1,6), (C 12 (I) ,1"  1 , 6) , (C44 (X) ,1-1,6), 

1  (RHO(l)  ,1=1,6) ,  (THICK  (I) ,1*1,6) 

51  FORMAT  (//',5X, 'Cl  1=  *  ,  6E  12.  3//,  5X,  *  C  12*  •,6E12.3// 

1,5X,'C44»  •,CE12.3//,5X,'BHC*  ',6112.3//,  5X , ' THICK-' ,6E 12. 3///) 
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c 

SBJECT 
sspuce  3 

C  CALCOLATE  A  PROCCCT  CP  LAIBt  IBASSPEB  HATtlCES, 

C  IN  PROD.  IN  WeAT  tCLlCWS.  SICRA  IS  TANGINTIAL  COnP  OP  INCIDEBT 
C  SLOWNESS;  ALPHA  If  PCRPEHDICUI A£  CCPPCklNT  CP  TBANSHITTEO 
C  LONGITUDINAl  SLCVNESS,  AND  EETA  IS  PERPENDICULAR  COHPONENT 
C  CP  TBANSniTTIL  SUIAB  SLCUHESS;  f,PP,Q,QP  ABt  LAIEB  THICKNESS 
C  IN  RADIANS  OP  ANGIE  -  SIN  AID  CCS. 

C 

DO  999  NNK  >  1.100 

THETA  =  EARSIK <(. ID-1) «PLOAI{NNN-1)) 

SICHA=K (1) *0 SIS (THETA) 

SIGSO=SIGPA**3 
DO  120  1=1, HF3 

120  ALPHA(I) =CDSQRT(  (K (I) •♦2-SIGSC) ♦ { 1 . ODO,  O.ODO)  ) 

DO  125  1=3, NF2 

125  BETA (I)  =  crS0RI ((KAPPA (1) •♦2-SIGSQ) • (1 . ODO.O- ODO)  ) 

DO  150  1  =  2, NP1 
TH=2.  *PI/K  (I)  ‘TfllCK  (I) 

P(I)=CDCOS  (A  1 1  HA  (r)*TH) 

PP(I)  =CCCCS(EIIA  (I)*TH) 

Q(I)  =  CDSIi;  (ALIKA  (I)  *TH) 

150  QP(I)=CDS1K(BI1A(I)*TH) 

CALL  CALCCP  (2,«AT1)  .  ' 

DO  151  1=1, NB1 
CALL  CALCCF  (1*2, HAT2) 

CALL  KATEDt  (BAT2,HAT1,PBOD) 

DO  152  11=1, « 

DO  152  JK=1,« 

152  aAT1(II,OP)=PECD (II.JK)  . 

151  COKTINUI 

DO  195  1=1,2 

DO  195  J=3,4  -  . 

PBOD(I,J)  =  FBCD (I,J) •OAI 

195  CONTINUE 

DO  196  1=3, «l  .  . . 

DO  196  J=1,2 

PBOD(I,J)  >  PECO  (I, J) *311 

196  CONTINUE  .  •  - 

C  . 

SEJECT 

SSPACE  3  ..  -  - . 

CALL  SIHEC (PNCD,ABP1,BBE1,ACPB)  .. 

XXB=ANP1  . 

XXI=ANP1*  (-JAT)  -  . 

TLn(NHH) =CIAES  (AHP1) 

TIP(NNK)=CATAK2(XXI,XXB) 

C  . . 

XXR=QNP1 
XXI=BNP1* (-JA1) 

TSH  (NNN) =CCAES (ENP1)  - . 

IP(TSfl(NNN)  .IT.  .IE-6)  CO  10  331 
TSP(KNN) =CATAX2(III,ZXB) 


331  XXB=AOPR 

XXI=AOpn*  (-JAI) 

BLH  (NNN) =CCAOS (AOPR) 
RLP(NNN) >LATAN2(IXI,XXR) 
999  CONTINUE 


113. 

11*1. 

115. 

116. 

117. 

118. 

119. 

120. 
121. 
122. 

123. 

124. 

125. 

126. 

127. 

128. 

129. 

130. 

131. 

132. 

133. 

134. 

135. 

136. 

137. 

138. 

139. 

140. 

141. 

142. 

143. 

144. 

145. 

146. 

147. 

148. 

149. 

150. 

151. 

152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 
161. 
162. 

163. 

164. 

165. 

166. 

167. 

168. 

169. 

170. 

171. 


88178(6. 9S0| 

950  rORtlKTCI  TBktSniTTBO  1CK6ITDDIMU  AB81IT00E  AlO  FBASB'/) 
DO  951  1-1.100 

88178(6.952)  71H  (I) ,  71t(I) 

952  POBHATC  •.3X.I12.5.3X.E12.5)  > 

951  COB7I80B  . .  ■ 

C 

8RI7E(6.953) 

953  POBHATC  1  TBASSHITTED  SBEAB  1BELI7DDE  AVD  PBASE'/) 

DO  955  1=  1.100 

Wfil7E(6,S5ll)  75n(I),  7SE  (I) 

954  POBHATC  •.3I.E12.5,3X.E12.5)  •  '  J  • 

955  COBTIHOE 

C  ■  .  * 

BRI7B(6.9S6)  "  - 

956  PORHATCI  BEPlECTtD  LO8CI70D1NAL  ABP1.I7D0E  ABD  PHASE*/) 

DO  958  1=1.100 

8RITE(6.9f7)  Rin(I)«  Etf(I)  . .  7- 

957  POBHATC  •.3X.E12.5.3X.E12.S) 

958  COBTIMUE 

C  .  ..  .  . . 

Stop 

EHO 


c  . . . .  . . .  ■■  ■  ■ 

SEJECT 

C 

C  SIHEQ  SETS  CP  A  SYSTEB  OP  4  ECOATIOKS  IB  4  DKKBCBHS 
C  ABO  USES  CBAHEB'S  BUIE  TO  SCIVE  FOB  3  OF  TBEB:  BEFLECTED 
C  AHPLITODE,  TBAKSHITTED  lOHCnODlNU  ABPLITOEE.  ABD 
C  TRASSHIITEE  SEEAE  AHPIITODE.  PHASE  IS  PBESENT  IH  THE 
C  COHPLEK  AHPtnODES.  THE  FOOETH  OHFKCWK.  THE  X  COHEOMEHT 
C  OF  PARTICLE  VilOCITI  IB  THE  SOLID  JCST  PAST  THE  LIQUID/ 

C  SOLID  INTERFACE.  IS  ACT  OF  IBIEBEST. 

C 

SOBROOTINE  SICEQ  (RATB , A AE.B AP. AOP) 

COnHON/CCAlC/llPHA(8) , EE1A (8) , P (8) ,PP (8) .Q (8) .QP(B). 
ICIKU)  ,C12(8),C4m8),K(e),KAPPA(8)  ,SIGHA,HP2 
COHPLEX*16  I!A1T(l»,4),HATf  («,4)  ,CC»ST(lJ)  ,TEHP  (4)  , 
1DETn,DET12,CET13,DET14,ANP,BNP.A0P.ALPBA.BETA, 
2P.PP,Q.CP.E 

BEAL*B  C11,C12.C44.X.I*AP.A.riGEA.ESlH 
C  -  . 


c 

c 


8SIH  =  2.*C44  (AP2)*SIGBA 
nATT(l.l)  =  SKHA/K(»P2) 

HATT(2.1)  =•  -  (C11(KP2)  »K(NP2)  -  2-*C44(8P2)* 

1  SIGnA*SIGEA/P  (BE2) ) 

BATT(1,2)  =  -  EETA{NP2)/KAEtA(»P2) 

flATT(2,2)  =  BSIR»HATT(1,2)  ,  . 

HATT(3.1)  =  A1EHA(AP2)/K1HP2) 

flATT(3,2)  =  SIGHA/KAPFA  {HP2)  . 

nATT(4,1)  =-  ESir*AIPHA  |HP2)/B(HP2) 

BATT(4,2)  »  •>C44{»P2)*(SIGHA*SIGHA-EETA(IIP2|*BET4(HP2))/ 
1KAPPA(NP2)  . 

HATT(1.3)  »  -FATB  (1,1)  >  .  . 

HATT(2,3)  *  -HATF(2,1)  * 

HATT(3,3)  =  -rATB(3,1) 

HATT{4,3)  -  -P.ATB(4,1) 


E  IS  THE  IBEECAACE  CP  THE  LIQUID; 
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172. 

173. 
171. 

175. 

176. 

177. 
176. 

179. 

180. 
181. 
182. 
183. 
18*1. 
165. 
186. 

187. 

188. 

189. 

190. 

191. 

192. 

193. 

194. 

195. 

196. 

197. 
196. 

199. 

200. 
201. 
202. 

203. 

204. 

205. 

206. 

207. 

208. 

209. 

210. 
211. 
212. 

213. 

214. 

215. 

216. 

217. 

218. 

219. 

220. 
221. 
222. 

223. 

224. 

225. 

226. 
227. 
220. 

229. 

230. 

231. 


c  r  IS  nvcEvi  or  ivcxoevcc  aicie. 

c 

E  «  C11(1)«K(1) 

r  ■  iLPHii  (D/r  (1) 

nETT(l(4)  *  K«Te(t.2)*e  *  n3TK(i,3)*p 

K(IT1(2,4)  =  HAIR  (2,2)  ♦  RATB(2,J)*r 

BATT(3.4)  »  CAIRO, 2)*B  ♦  8A1B(3,3)*P 

BATT(4,4)  nATR(4,2)*E  ♦  nA1B(4,3)*r 

C 

C  THE  BAT8IX  Cf  CClFflCIEHIS  OF  THE  SET  OP  4 
C  EQUATIONS  IN  TFE  4  UBFNCHNS  ANF,  ENF,  AOP  IS 
C  HON  BEAD!.  NEXT  PUT  TEE  COEETAN1S  IN  THESE 
C  EQUATIONS  IN  4-1ECTCB  CONST  (I). 

C 

COBST(I)  =  -EATR  (1,2J •E»HATB(1,3)  *2 
COHST(2)  =  -KATR  (2,2) •E«nATf (2,3) »P 
CONST(3)  =  -flATR(3,2)  •E*flA7B(3,3)  •? 

COUSTIU)  =  -EA18(4,2) •E*BA1E(4,3) 

C. 

C  VON  USB  CRAHEB'S  BUIE  TO  SOLVE. 

C 

CAtl  COET4  (BATT.DETII) 

DO  731  1  =  1,4 
TEHP(I)  =  BATl  (1,1) 

HATT(I,1)  =  CCkST(l) 

731  CONTINUE 

CALL  CDET4  (HATT,CET12J 
AHP  =  0ET12/CET11 
DO  732  I  =  1,4 
HATT(I,1)  *  IKBPd) 

TESP(I)  =  HATT)I,2) 

BATT{I,2)  =  CONST  (I) 

732  CONTINUE 

CALL  CDET4(RAn,0ET13) 

BNP  =  DET13/tIT11 
DO  733  I  =  1,4 
HATT(I,2)  =  TII!P(I) 

TEHP(I)  *  BATT(I,4) 

HATT(I,4)  -  CCKST{I) 

733  CONTINUE 

CALL  CDET4(HATT,riI14) 

AOP  =  DET14/DET11 
BETOBH 
BHD 
C 
C 

SEJECT 

C  CALCOF  CALCCIATFS  THE  TBANSBEB  BATBIX  OF  TBB  ZTH  LA7EB 
C  AND  PUTS  THE  BESOET  IN  A.  SEE  IBBONS*  EQ.  5.10 
C 

SUBROUTINE  CAICOF  (I, A) 

COflflON/CCAlC/)lPHA(8)  ,EETA (6)  ,P(e) ,PP (B) ,Q (8)  ,QP(8)  , 
1C11  (8)  ,C12  (8)  ,C4  4  (8)  ,K  (B)  ,i(APPA(8)  ,SlCnA,NP2 
BEAL*8  C11,C12,C44,SICnA,K,EAFEA,CA,CE,CC,SICSQ,KA,KSQ, 
lKAr,KAPSC,tirFC,CKSQ,CKAfSC 

COBPLEX*  U  AIFNA,CETA,F,FF,C,CC*A (4,4) , ALPH, AIPHSC* BET, 
IBETSQ, PA,FFA,CA,CPA, DlEfSC,CSCASQ 
CA>C11(I) 

CB>C12(1) 

CC«C44 (I) 
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233. 

23%. 

23$. 

236. 

237. 

238. 

239. 

240. 

241. 

242. 

243. 

244. 

245. 

246. 

247. 

248. 

249. 

250. 

251. 

252. 

253. 

254. 

255. 

256. 

257. 

258. 
259*. 
260. 
261. 
262. 

263. 

264. 

265. 

266. 

267. 

268. 

269. 

270. 

271. 

272. 

273. 

274. 

275. 

276. 

277. 

278. 

279. 

280. 
281. 
282. 

283. 

284. 

285. 

286. 

287. 

288. 

289. 

290. 

291. 


siGso-sicnk**3 

ALPH^ALPHJ)  (I) 

&LPHSQ>AL(I:**2 
BEli^BBTA  (I) 

B2TSQ-BET*«3 

ICA>K(I) 

KSQi>KA**2 

KAP>KAPFA(I) 

KAPSQ*KAP«»2 

DIPPC-CA-CB 

OIfrSQ=SICSC-fEISC 

PA=P(I) 

•  PPA-PP(I) 

QA=Q(I| 

0PA=QPCX) 

CKSQ=CA*KSQ 

CKAPSQ=CC*KAPSC 

CSCASQ=CB»£IGJC*CA*A1PBSC 

A(l,1)=  (C1F1C«£1GSC*PA»CSCASQ*PFA)/CKSQ 

A  (1,2)=SlGnA/CKSC* (EPA-tA) 

A (1,3) =-SIC«A*CIFFSQ/(ArFB*KAFSC) ♦C A-2. •BET*StG«A/XAPSQ*QPA 
A  (1,a)=-SICSC/ (ALPH*CKAPSQ)*0A-EET/CKAPSQ*QPA 
A  (2,1)=CSCAS0*tlFFC*SIGeA/CKS0*(PPA-PA) 

A (2,2)  =(C£CASC*PA*DIFFC*SIGSO*PFA)/CKSC 

A  (2,3) =CSCASC*tIFFSQ/ (ALPH*KAPSC) •0A-2*CIFFC*SIGS0*BET/KAPSC*QPA 
A  (2,4)=CSCASQ»SI GH A/ { ALF K*CKA PSC)  •CA-DItFC*SIGHA*BET/CKAES(.*CPA 
A(3,1)  =DtFFC*AtPH*SIGBA/CKSC*CA-CSCASC*SIGflA/(B£T*CKSe) •QPA 
A (3,2) *-AlEH/CKSC*CA-SIGSQ/ JEET^CKSQ) ‘QPA 
A  (3,3)=  (-D1FFSC*FA*2.*SIGSC»PPA) /KAPSQ 
A(3,4)  *SIGKA/CFAPSC*(FPA-PI) 

A  (4,  1)=-2.  ♦CC*riFFC*AlPl)*SJGS0/CKS0*0A*CC*CSCASQ*DIPPSQ/(BBl*C3tSQ) 
1*QPA 

A (4,2) *2.*CC*AIPn*Sia«A/CKSC*CA*CC*CIFPSC»SICaA/(EET*CKSQ)4QPA 
A  (4,3)=2.  «CC»riFFSC»SIGF,A/KAPSQ*  (PA-PPA) 

A  (4,4)*  (2. •SIGSQ  =  PA-D1FFSC*PPA) /KAPSQ 
RBTUBH 

•  ERO  ■ 

C 

SEJECT 

C  BATnOL  BOLTIPIESEATRZCES  A  AtC  B,  PDTTIB6  2BB 
C  BESULT  IB  C. 

C 

SOBBOOTIHE  HAZeUt  (A,B,C) 

C0flPLEX*1£  C  (4,4),A(4,4),B(4,4) 

DO  100  1*3,4 

DO  100  0=1,2  • 

A(1,J)*-A(I,J) 

100  B(I,J)  *-B  (1,0) 

DO  150  1*1,4  . 

DO  150  0*1,4 
C(I,J)*(0. 0,0.0) 

DO  150  K=1,4 

150  C(I,J)=C  (1,0)  «A(I,K)*B.(K,0} 

DO  200  1*3,4 

DO  200  0*1,2  . 

A(I,0)— A(I,0) 

B(I,0)*-E  (1,0)  » 

200  C(I,0)  =>C(I,0) 

BETOBB 

END 


I 


292. 

C 

CAIDET  CALCDLITtS  THE  OETEBelRAII  01  A  3l3  BATBIf  D. 

293. 

C 

PUTTING  THE  filSUll  IH  OET. 

2S9. 

c 

29S. 

SUBROUTINE  CAIEET  (D.DEI) 

296. 

C0nPLBI*16  C(3,3) ,OET 

297. 

DET<’0(1.  1)  *D  (2,2)»D(3,3»»D|I,2)*D(2,3)*D|3, 

1)  *0(3,2)  • 

298. 

1-0(3,  1)  *C  (2,2)  *0  (1,3)-C(2,  1)*0  (1,  2)  *0  (3,3)- 

D(3,2)* 

0(2,3)*D(1 

299. 

RETUEE 

aoo. 

EEO 

301. 

c 

.  •  • 

.. 

302. 

SSP4CE  3 

303. 

C 

304. 

c 

C0ET4  CAtCOIATES  A  4X4  DETEBBIHABT  BX  EXPAMSIOI  BT 

305. 

c 

BINORS  ON  THE  fIRST  CCIOaN. 

306. 

c 

.  . 

307. 

SOBROOTINE  CDET4(D4,DET) 

308. 

C0RELBX«1£  C4(4,4),  0£1.  03(3,3),  TOET 

» ■ 

309. 

c 

• 

310. 

DO  831  1=1,3 

- 

311. 

DO  831  3=1,3 

312. 

D3  (I,J)  =  C4(14l,J«1) 

313. 

831 

CONTINUE 

314. 

CALL  CAICIT (t3,TCIT) 

315. 

DET=TDET*D4(1,1) 

316. 

DO  832  J=1,3  . . 

317. 

03 (1,J)  =  C4  (1,0*1) 

318. 

832 

CONTINUE 

319. 

CALL  CALDET(t3,TCET) 

320. 

DET  =  OET  -  TCET*D4(2,I) 

At. 

321. 

DO  833  0=1,3 

322. 

D3(2,0)  =  04(2,0*1)  -  t  .  . 

-  -  •  ■  • 

323. 

833 

CONTINUE  i 

.. 

324. 

CALL  CAI0ET(D3,TDET) 

325. 

DET  =  DET  ♦  TtET*D4(3,1)  . . . 

» 

,  . 

326. 

DO  834  0=1,3 

t 

327. 

D3(3,0)  =  04(3,0*1) 

Jt- 

328. 

834 

CONTINUE  . . 

'f,:.  ’ 

*• 

329. 

CALL  CALIET  (C3,TCET) 

330. 

OET  =  DET  -  TCET*04X4,1) 

. 

331. 

RETUBH  . -  -  . 

332. 

C 

0- 

333. 

C 

CLASS  DATA: 

334. 

c 

0.5760CD12  0.20200D12  0.21600D1 

0.25D0 

335. 

c 

PEREX  DATA  (AUIC) : 

336. 

c 

7.30E11  2.50E11  .232B1 

• 

1B1 

337. 

c 

HEAVE  PEREX  EATA 

338. 

c 

7.30E11  2.50E11  .280E1 

.1E1 

339. 

c 

HEAVr  SILICATE  ELINT  GLASS  CATA  (ADLD) : 

340. 

c 

6. 13E11  2.18E11  .3879E1 

4» 

1B1 

341. 

c 

T-40  GLASS  CATA  (AULO) : 

342. 

c 

6.30E11  2.26111  .339B1 

• 

181 

343. 

c 

LIGHT  BORATE  CBCNH  GLASS  DATA  (ADLD): 

344. 

c 

5.82E11  l.aiEII  .2243E1 

.181 

345. 

c 

HATES  DATA: 

346. 

c 

2.2SE10  O.OOEOO  1.00E00 

1.00E00 

347. 

c 

ALUnlNUn  CATA: 

348. 

c 

1. 11D12  2.50011  .2695B1 

• 

3SS0 

349. 

c 

SAPPKIRC(ASSOnTC  ISOTROPIC)  CATA: 

350. 

c 

4.94C12  1.45C12  3.98C00 

I.OOEOO 

351. 

c 

CHROBIUn  DATA: 

-  83  - 


3S2. 

C 

3.50E12 

1.01E12 

7.20C00 

378-1 

353. 

c 

POlt  SI  0411: 

351. 

c 

1.66012 

.800012 

.233E1 

355. 

c 

CU4RTZ  Dill: 

356. 

c 

.87012 

.58012 

.26501 

.04X00 

356.1 

c 

30SZD  S11K4  C4T1: 

356.2 

c 

7.85011 

3.12011 

.22001 

356.3 

c 

COPPER  Dill: 

356. « 

c 

16.9011 

7.54011 

.8901 

356.5 

c 

COPPER/9.98  Its 

41 

356.6 

c 

16.0011 

7.66011 

.75801 

356.7 

c 

RLORIHOn/IC  S  CO 

lllOr  122  ClTl 

i: 

356.8 

c 

11.7011 

3.0011 

.29501 

.3580 

357. 

c 

POLYCTHYLEKf  0174 

J 

358. 

c 

3.40E10 

.260E10 

.900E00 

1.080 

359. 

c 

•Y4Ca0H'  0474: 

360. 

c 

1.0  E2 

1.0E0 

1.0E-6 

.18-3 

361. 

c 

flTRIRITB  SC  Elf 

*  83.5  (CILC. 

EEBKOlilTZ)  Dill 

• 

362. 

c. 

5.7E10 

t.43E10 

1.3 

*1.080 

363. 

c 

TITBIMITE  SC  Clf 

>  88.8  (C4LC. 

BEBKOklTZ)  04T1 

• 

• 

364. 

c 

7.73E10 

1.S3E10 

t.3£C0 

1.080 

365. 

c 

EXIBITE  SC  OAF  » 

83.5  (ClIC.  BEBK0EI7Z)  01T4: 

366. 

c 

3.27E10 

8. 18X9 

1.2E0 

1.080 

367. 

c 

EXIBITE  SC  C4P  * 

88.8  (ClIC.  BE1P0HI7Z)  C1T4: 

368. 

c 

8.  2  IE  10 

2.C5I10 

1.3B0 

1.080 

369. 

c 

CIRBOB  V»3  IlTl: 

(2*7) 

370. 

c 

2.1E11 

5.24C10 

.23311 

1.080 

371. 

c 

CIRBOB  V»4  ClTl: 

(2-7) 

372. 

c 

2.  BE  1 1 

7.CE10 

. 175E1 

1.080 

373. 

c 

CIRBOB  V»5  ClTl: 

(2=7) 

374. 

c 

3.5E11 

8.7SE10 

.140E1 

1.080 

375. 

c 

C4RB0H  Z«6  ClTl: 

(7.4) 

376. 

c 

2.4E11 

6.0E10 

.1581 

1.080 

377. 

c 

CIRBOB  Zs'S  ClTl: 

(»*«) 

378. 

c 

3.2E11 

8.0E10 

.20B1 

1.080 

379. 

c 

CIRBOB  Z-10  0171: 

380. 

c 

4.0E11 

1.0B11 

.25E1 

1.080 

381. 

c 

PLlTIBOa  ClTl: 

381.1 

c 

2.27012 

6.4011 

2.1401 

.88-1 

381.2 

c 

COLO  04Tl: 

381.3 

c 

2.07012 

2.85011 

1.9301 

.15980 

381.4 

c 

TlTlBIun  ClTl: 

381.5 

c 

1.66012 

4.4011 

4.5000 

.210-1 

382. 

EBP 

383. 

$0614 

384. 

2 

386. 

2.2SE10 

O.OOEOO 

1.00E00 

1.00800 

388. 

5.B2E11 

i.eiEii 

.2243E1 

.7780 

390. 

7.30E11 

2.5CE11 

,232E1 

.081 

392. 

7.30E11 

2.50C11 

.232E1 

.181 

400. 

SSTOP 

> 

I 


0.1 

0.2 

O.J 

0.« 

0.5 

0.6 

0.7 

o.e 

0.9 

1. 

2. 

3. 

«. 

5. 

6. 
7. 
6. 
9. 

10. 

11. 

12. 

13. 

1«. 

15. 

15.1 

16. 

16.1 

17. 

18. 
19. 

19.1 
21. 
22. 

22.1 


32PMD12  8 

PSOCRAH  TC  CAICOIATE  THE  V(Z}  RESPONSE  OF  AN  ACOUSTIC 
LEVS  FOB  A  SAHFIE  WHOSE  REFIECTAVCE  IS  CITEV  AS  DATA. 


//TZ  JOB  •F17]E7,181S'B  EBA1*,CIASS«Z,BE6ICN>S12X 
/•JOBPAfifl  DEST=SI1F 
//  EXEC  WATFIV 
//GO.SYSIH  DD  • 

SBATFlV 

C  THIS  PROGBAf  CAICULATES  TEE  T(2)  RESPONSE  CP  AH  ACOUSTIC 
LENS  TO  A  TARCET  KKOSE  REFLECTANCE  FUNCTION  IS  SPECIFIED 
AS  DATA.  THE  PUPII  FUNCTION  CF  THE  TENS  AND  THE  BACK 
POCAL  PLANE  CISTRIEOTION  MUST  ALSO  FE  GIVEN  AS  DATA. 

OTHER  INPUTS  ARE  ACOUSTIC  FRECUENCV,  LENS  BACIUS  AND 
OPEUING  ANGLE,  2  EINIKUH  ANC  BAXIHUM.  NATER  IS 
ASSUMED  AS  CCUPLING  HEDIUH  WITH  TERIERATUDE  60C. 

THIS  PSOGRAE  NUnEElCALLI  EVALUATES  THE  V(Z)  INTEGRAL 
IN  A.ATALAR'S  THESIS  <STANFCRD  UNI VERSITY ,  1978)  GIVEN 
AS  EQ.  (U.26).  IN  ADDITION.  THE  VAEIATICN  IN  PATH 
LENGTH  WITH  ANGIE  FRCfl  LENS  TO  FOCUS  IS  INCORPORATED. 


REAL*4  ZG  (101) ,VZ10G{101) 

REALMS  PI,F8EC,RIENS,OFANG,2tlIN,ZtlAZ,RLa,5LP, 
1U1X,01Y,U1A,U1P,PA(81) ,PF(81) ,E AP, F ,ZI NC,KO, 
2AREA,R,VZA,V2X,V2Y,VZP.ANG,EX,FY,PX,PY,Z,ALPBA, 
3Q1,Q2,RB,e£lAR,ANG1 
COMPLEX* 16  JAY,BEFL(100),P(81) ,E?Z, 
1XEXP,VZ(101) ,U 1(201)  ,VZHAX 


23. 

PI=3.1415S27 

24. 

JAY« (0.0C0,1.0C0) 

25. 

C 

26. 

C  PBEQ  IS  ASSOHED  TC  BE  IN 

GHZ; 

27. 

C  IN 

HICPONS;  OPANG  IN  DEGREES; 

28. 

C  INCREHENTS (.01)  Cl  SIN(THETA) 

29. 

C  THE  LENS  APERTURE  BACIUS; 

f  I 

30. 

RBAD,FBEQ 

31. 

READ,  Rl£KS,CFANG 

32. 

READ, zniN, ZBAX 

33. 

HRITE(6,222)  FREQ 

34. 

222 

FORKAT(»1  FFEQUENCr 

s  f 

35. 

HRITE(6,223)  RIENS, OPANG 

36. 

223 

rORHATCO  LENS  RADIUS  * 

37. 

1,  E10.4,  •  DECREES*) 

38. 

HRITE(6,224)  ZRIN,  ZHAZ 

39. 

224 

FORHATCO  ZriN  ■  *, 

E12 

40. 

1*  niCRONS*) 

41. 

DO  10  1^1,100 

42. 

READ, Rin, RIF 

43. 

RX>=RLn*DCCS  (RIF*  1.  DO) 

44. 

RY=BLH*0SIN(Brr*  1.00) 

45. 

REFL(l)»[CEriX(nx,RX) 

45.1 

10 

CONTINHE 

46. 

HRITE(6,22S)  REP1(1), 

REP 

47. 

225 

rOKHATCO  CifL(l)  - 

BLENS  AND  ZHIH 
REFLECTANCE  IN 
01 (100)  IS  01 


AND  ZHAZ 

EQUAL 

AT 


1  DEGREE  INCBEHEHTS. 

E10.H,  •  CHZ») 

',E10.4,'  RICB0N5;  OPENING  ANGLE 


a  t 


ZBAX  -  Bl2.4t 


.5,  IX,  B12.5, 


-  85  - 


REPL(IOO)  -  *, 


I 


«8. 

1  812.5.  1Z.  £12.5) 

DO  20  i«1.201 

SO. 

51. 

■EAD.UIA.UIP 

52. 

Oil  «  D1A«I)StK(aiE*1.D0) 

53. 

018  >  01A*tC0S(U1P*1.00) 

5«. 

01(1)  -  DCEPtX (U1X.01Y) 

55. 

20  COHTIHOE 

5«. 

VR1TS(6.22£)  01(1).  01(201) 

.57. 

226  FOSflATCO  01(1)  *  E12.5. 

M,  E12«5«  •  01(201) 

58. 

1  812.5.  IX.  £12.5) 

59. 

BEAD(S. 16)  (PA  (X) .K=1.81) 

60. 

B8A0(S,16)  (PP(I<)  ,K=>t.81} 

61. 

16  POBflAT<1H,5ri2.8) 

62. 

DO  21  I»1.81 

.. 

63. 

PX«PA(I)  ♦CCOS  (PP(I)*1.D0) 

68. 

PTsPA(I)  •ESIN  (PP(I)  *1.00) 

. 

65. 

P(I)=DCKPIX(PX,PY) 

66. 

21  CONTIUUE 

.. 

67. 

«RITE(6,227)  P(1),  P(81) 

■ 

68. 

227  P08nAT(<0  P(1)  «  •,  £12.5.  1Z,  812.5.  •  8(81)  > 

69. 

1  E12.5.  IX.  812.5) 

• 

70. 

C 

- 

" 

71. 

c 

72. 

BAP=BLERS  «  CS I E (0PAHG*PZ/180.) 

73. 

P«1.15*RLE1IS 

* 

78. 

II1IC={ZBAX-2H1K)/100. 

75. 

K0=2*PI*FfEC/1.5 

76. 

c 

• 

77. 

c 

1.5  niCRONS  PER  RANCSECORD.  78AX 

IS:  SPEED  XB  BATES 

78. 

c 

CHARGE  SPEED  ECB  AHCtHEB  IICUID. 

IIKEHISE 

79. 

c 

.0102  NSEC*NSEC/KICRCK  IS  A1TE80AT20N  COBSIABT  IN  BATES 

79.1 

c 

AT  60C  TEHPERATCBE.  CHARGE  EOS  AKCTHEE  IICUID  OB 

79.2 

c 

AHOTHIR  TEHFEBATUBE. 

80. 

c 

....  ...... 

80.1 

ALPHA«.0102*FEEC*PREQ 

. 

80.2 

c 

. 

81. 

8BJECT 

... 

82. 

c 

BHAT  FOLLOHS  IS  A  CCOOLE  LOOP  TO 

CALCOLATE  V(Z). 

83. 

c 

THE  IHBER  LOOE  (CO  100  1=1.101)  CALCULA1ES  1Z(J)  FOR 

88. 

c 

A  FIXED  VALUE  CF  Z  HIIICU  IS  IRCSEBERTED  IB  THE  OOTEB 

85. 

c 

LOOP  (DO  120  0=1.101). 

86. 

c 

87. 

Z>ZHIII 

87.1 

yZBAX  «  (0.00,0. DO) 

88. 

DO  120  0=1,101 

89. 

z=zniv«DrioAT  (o-1)*zibc 

... 

90. 

ZG(0)  =  5RG1(Z) 

91. 

c 

92. 

c 

Z6(J)  IS  SXBPIY  A  TABLE  CF  Z  FOB 

ELCtlllG  BOOTHES 

93. 

c 

• 

98. 

TZ(J)  >  (O.OO.O.DO) 

9 

95. 

c 

.  t 

96. 

00  100  1=1,101 

f- 

97. 

B-  OFLOAT  (I>1)/100. *882 

98. 

AREA  «  2«PI*R*8AP/100. 

*  *  .  .. 

99. 

DVZ-01(I)*U1(1)*ABZA 

. 

100. 

c 

. 

101. 

c 

P  CORES  AS  A  FOMCYICM  OF  ARCIE  XB 

EEGBIES.  CONfEBT  TO 

102. 

c 

THE  APPBOPBITE  VAIOE  CP  RADIOS. 

. 

102.1 

c 

MOTE  THIS  VERSJCM  EXPECTS  A  tOPXX  FDBCTXOE  18 

86 


t02.2  C  t  OeSREE  tKCEfRENTS,  VHEREES  CtEER  fBRSIOMS  EXPECT 

102.3  C  2  DEGREE  INCSERENIS.  OU1PU1  fSOH  PUPIL. EO  IS  SUE  TABLE. 

103.  C 

10«.  AKC1>=DARSlN(R/(RtENS)) 

105.  XAMG1«lDlkl (A SCI*  ISO. /3. 141 59) ♦I 

106.  DVZ-DVZ*P  <UKC1)  *P(IAHG1) 

107.  C 

108.  C  BEPLECTANCE  CCHES  AS  A  POKCTICX  OP  SIK  OP  INCIDE8CB 

109.  C  ANGLE.  PIMO  AIIFCFRIATE  VALUE  fOR  A  GIVEN  RADIUS. 

110.  c 

111.  XAV6=IDINI  (<R/(1.  15*BLEliS))*100.)  »1 

112.  DVZ-DVZ«REFL(XANG) 

113.  C 

114.  XEXP=CDEXF (JAX*2.*(K0*JAI*4tPHA)*Z* 

114.1  1DSQHT(1.- (R/f)*(B/P))) 

115.  DVZ°DVZ*ZEXF 

11S.00S  C 

115.01  C  THE  NEXT  LINES  ACCOUNT  POP  VAEIAGLE  VATEB  ATTBBDATIOI 
115.02  C  OH  THE  WAX  IC  THE  PCCOS 
115.03  C 

115.04  BR  «  RLENS*0CCS(AKG1) 

115.05  BSTAB  s  tSCRT ( (. 15*RLERS*BE) • (. 15*RLENS*RB)  »  B*B) 

115.06  ANG  =  DAGSIN (E/BSTAR) 

115.1  DVZ-07Z*D£XE(AlFIIA*2.*t1.15*BLENS-ESTAB)) 

115.11  C 

116.  VZ(J)=V2(0)*DV2 

117.  100  CONTINUE 

117.1  C 

117.2  C  HON  KEEP  TRACK  Cf  LARGEST  HAGHIIUEE  TO  NOfinALZZS  LATER. 

117.3  C 

117.4  Q1  *  CDAES(VZ(J)) 

117.5  02  «=  CDABS(V2I1AX) 

117.6  ZP(Q1  .GT.  02)  VZHAX  *  VZfJ) 

118.  120  CONTINUE 

119.  C 

119.1  C  PRINT  OUT  NUEEBICAt  VALUE  CP  TZflAX 

119.2  C 

119.3  NRITE(6,2€9)  VZHAX 

119.4  269  FORHATCO  VZHAX  -  E12.S,  U,  £12.5) 

119.5  02  *  CDABSIVZEAX) 

119.6  VRITE(6.271}  C3 

119.7  271  rORHATCO  HAGNITUDE  OE  VZHAX  •  E12.5) 

120.  SEJECT 

120.1  C 

120.2  C  NON  NORHALIZE  RELATIVE  TC  BAXIHUB  RETURN. 

120.3  C 

120.4  c****************** 

120.5  C  DEFINE  VZHAX  AS  GLASS  VZHAX  TC  COHPARE  IIFEERENT 

120.6  C  VZ  CURVES  CK  SAHE  SCALE.  Tb'IS  VZHAX  IS  FOB  A 
120.65  C  24  HICROU,  66  CEGREE  lENS  C  'LAYER*  TARGET. 

120.7  c*************************************************** 

120.9  VZHAX=<0. 1588CD-01,0.543e0D-01) 

121.  HRITS(6,1S01  r 

122.  190  FOBHATCI  V  (Z)  ARPIITUDE  IN  DB  AND  PHASE') 

123.  DO  200  I»1,101 

123.1  VZd)  «  V2(I)/VZaAX 

124.  VZA=20. *010010  (CCABS(VZ(X))> 

125.  TZLOG(Z)  >  SNGl(VZA) 

12b.  VZX>VZ(I) 

127.  VZy«VZ(I)*(-JAt) 
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128. 

TZP*06TkN2(VZI.TZT) 

129. 

«RI1'E(b,201)  VZA.VZP 

130. 

201  rOHHATC  •«2X.et2.S.3Z,E12.5) 

131. 

200  COKTIHUe 

132. 

c 

133. 

C  CAIL  PLOTTZtC  BCOTIMES 

134. 

C 

139. 

CALL  STABTG('CIMIL.LliiE8Z«e**,0 

.0) 

135.1 

.  CALL  StlBJIG<-IO..-45.,£.,IS.) 

135.2 

CALL  SETSfGC 132,1.) 

136. 

CALL  GDAPtCCVEBIT**,  1.101, ZC.VZLOG.O.'Z*',' V(Z)  IM  DB* 

137. 

1'VITfiIHllE  ZDAF  -  83.5 

V(Z)  FS2.6  GHZ*') 

138. 

CALL  EXITG 

139. 

STOP 

140. 

BIO 

141. 

SOATA 

142. 

2.6 

143. 

24.  66. 

144. 

-10.  5. 

145. 

0.289460  00  -C.59930D-17 

146. 

0.269440  CO  0.899020-17 

147. 

(  IR  ALL,  100  TAtOES  OF  REFLECTANCE  AHPLITODE  C  PHASE.] 

148. 

r  THIS  IS  OC1PU1  CF  REFLEC.SIA  EACGBAB.  ] 

243. 

0.603300  00  -0.313650  01 

244. 

0.698710  00  -0.314060  01 

245. 

0.23119E-02  -0-35547E  01 

246. 

0.23117E-02  -0.3S553E  01 

247. 

[  IN  ALL,  201  VALUES  CF  01PI0S.  E.C. 

FBCB  AlAIAR'S  PBOGBAH. ] 

444. 

0.86933E-03  -0.39S22E  01 

445. 

0.66643E-03  -0.39SS6E  01 

446. 

3.5592  3.5597  3.5612 

3.5637  3.5672 

447. 

3.5718  3.5775  3.5642 

3.5921  3.6012 

446. 

r  IN  ALL,  81  VAIOES  OF  POFIt  FOICIICIi  ABEIITUDE, 

449. 

E.G.  LENS  TBANSniSSION  FUNCTION 

FBOB  LEBONS* 

450. 

PBOCRAH,  RCDIFIED  AS  P0EIL.80. ] 

461. 

0.7489  0.7084  0.6686 

0.6291  0.5897 

462. 

0.5499 

463. 

-0.6708  -0.6713  -0.6728 

-0.6753  -0.6788 

464. 

-0.6833  -0.6869  -0.6955 

-0.7031  -0.7118 

465. 

(  IN  ALL,  81  POFII  FUNCTION  PHASES 

TO  GO  KITH  A0PLITODE. ] 

478. 

0.1830  C.1764  0.1701 

0.1640  0.1578 

479. 

0.1515 

500. 

ZSTOP 
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IPPZIOIX  c 

PBOCBAH  TO  CALCOLATE  THE  PROTCACCDSTZC 
BBSPOMSE  III  A  I-BEGIOM  SISTER. 


//OPTAC2  JOB  'FlISEIflBI'.'B  EBAI',CI.ASS>E,BECIOa>S12K 
/MOEPABB  DEST=SEtP 
//  EXEC  KATPIT  kOLlST 
//GO.SYSIM  DD  • 

SBATPIV  NOLIST 

C  THIS  PORTBAN  PBCGRAH  SOLVES  FOB  THE  ACOCSTIC  POVEB 
C  RADIATED  INTO  TIO  SEEI-INFIKITE  HAIP-SFACES,  BEGIOHS 
C  1  AND  2,  AND  GENERATED  BY  OPTICAL  ABSOBBTION  IN  BEDIUB 
C  2.  THE 

C  INCOniNG  OPTICAL  POHEB  IS  ASSOBED  TO  BE  INCIDENT 
C  THROUGH  TRANSPABEBT  BEDIUB  1,  AND  IS  ASSURED  TO  BATS 
C  UNIT  INTENSITY, 

C  THE  ACOUSTIC  GENERATION  IS  SCtVED  BY  THE 
C  HETHOD  OP  COOPIEI  fiCCIS  IN  A  DISTRIBUTED  FINITE 
C  SOOBCE  REGION  AS  PRESENTED  IN  £. A.  AOIC'S  TEXT  V.I,  CB. 6. 
C  THE  PRESENT  ERCGRAB  MAS  HBITTEN  BY  BC  BEAT  IN 
C  SEPTEBBEB  1960. 

C 

DIHENSION  RH0(2),  C (2)  ,  C11{2),  BETATB(2),  AtPHA(2), 
1A(2),  Z(2),  TAUL(200),  STR  (200) 

BEALK(2)«  XE(2),  B,G, R0,B1,L,LKAX,LINCB,S1 (20  1) , 

1  S2(201) ,F1(2C1) ,F2(201) ,PB1N,FEAX,E(201) 

COBPLEX  TAU(2),ZTn(2) ,  R,TE1 1 , TE20 , TH22, 
1JAX,UU,VV,££TA,AF2F,A«1F,AP20, ABIO 

c 

C  PBEQ  FRECUENCr  OF  BODDLATION  OF  IICUX 

C  BET  «  OPTICAL  ABSCRETICH  COEFFICIENT  OF  2 

C 

BEAD,  FBIN,  FBAX 
HRITE(6,133)  FtllR,  FflAX 

133  rORBATCI  FHIA  =  •  ,E1 8.8,3Z,' FB AX  -  E18.8//) 

READ, BET 
VRITE(6,50)  BET 

50  FORBATCO  OPTICAL  ABS  COEFf  BBT  »  E20.8,*  1/H') 

BETA  -  CEPIX  (BET,0.) 

HRITB(6,3C) 

30  FORBATCO  KK  RBO 

1  C11  BBTATH*) 

DO  45  I’>1,2 


XX(I)  B  THEFBAL  COBCOCTIVITI  OP  ITB  BECIOH 

iiHO(I)  -  DENSITY 

C(I)  >  SPECIFIC  HEAT 

011(1)*=  ELASTIC  STIFFNESS  ClI 

BETATR(l)  >  THFRHAL  COEFFICIENT  OF  LIHEAB  EXPABSIOI 

READ,  KK(I),  RHO(I),  C(I|,  C11(I),  EETATR(I) 

HRITE(6,4C)  KP(I),  BHO(I},  C(l),  C11(l),  DETATH(I) 

)  FORnATCO',SE2C.8) 

)  CONTINUE 

SET  OP  LOOP  TO  CALCULATE  BESOITS  FCR  201  VAIOES  CF  FBEQOBNCr 


PWI.  1 11 IJ...  II  k  M  ■  I 


70.2 

70.3 

70.4 

70.5 
70.7 

73. 

74. 

75. 

76. 

77. 

77.1 

77.2 

78. 

79. 

80. 
81. 
82. 
63. 

84. 

85. 

86. 

91. 

92. 

92.5 

93. 

98. 

99. 
106. 

107. 

108. 

109. 

110. 
111. 
112. 

113. 

114. 

115. 

116. 

117. 

118. 

118.5 
119. 

120.5 
121. 

121.5 
122. 

122.5 

123. 

123.5 

124. 

124.5 

125. 

125.5 

126. 

126.5 

127. 

127.5 

128. 

128.5 


C 

riNCB  -  (61CG1C(FH4K)-4tOC10(rniH))/200. 


DO  73  H 

»  1,2C1 

FBBQ  - 

ALCCIO(fniB)  *  (H-1)*FIHCB 

FBEO  - 

10.**FBEQ 

• 

F(H)-FREO 

c 

CALCULATE 

ALPHA  (1) 

-  TKEBRAL  DIFFDSIVITY 

or  ITU  LiTEI 

c 

A  (I) 

«  TUEBEAL  CIFFDSIOB  CCEFFICIEBT 

c 

TA(((I) 

=  (1*J)*A(1) 

c 

Z(l) 

-  ACOUSTIC  IHPEtASCE 

c 

K(I) 

-  ACOUSTIC  inOPACATlOH 

CO AST A NT 

C 

ZTB(l) 

«  THBBHAL  inPECANCE 

c 

JAY-(0. 

»1-) 

DO  60  I 

*1,2 

ALPHA(I)  =  KK|I)/(BHO|I)*CjI)) 

A(I)  =  SQB1((2.*3.14159*PHEC) /<2.*AtPHA|I))) 

z{i)  =  saBT(niic(i)*cn  (D) 

T*0(I)  =  CKPIX  lAdt.Ad)) 

K(r|  =2.*3. 14159*FBBC*SCB7(8HC{I)/C11(I)) 

2TH(I)» (l.*JAY) *SC8T(3.14159*IfilQ*BBO  (I)*C(I) 

1*KK(I)) 

60  COUTIHOE 

TSD  *=  l./A(2) 

C  VHITE(6,«*)  TSD 

C  65  FOBHATCO  TEIRBAL  SXIB  DEETB  A(2)  =  B20.8) 

C 

CAli  TBAPS(1, 1,10000, 1,1) 

C 

C  CALCOLATE  TRII,  TH20,  AMC  T1I22,  lERPEBATUBE 
C  ASPITUOeS  AX  IBE  Z«0  BODiilDABI.  .  4 

C 

TK11  *=  BETA/(2TB  (1)4(1. ♦Z1H(2J/ETH(1))4(BETA*TA0{2))) 
C 

TH20  =  BETA*  (BETA* (ZTH(2)/Z1H(1))*TA0{2)) 

TH20  =  1H2C/(2TH(2)) 

TH20  =  1B20/(1.*ZTH(2)/ZTH(1)) 

TH20  =  TH20/ (EETA*BETA  -  TAU(2) «TAD (2) ) 

C 

TH22  •=  BE7A*TA0(2)/ZTH  (2)  - 

TH22  »  TH22/(EEXA4B£TA-1A0|2)4XA0(2)) 

IH22  -  -TJi22 

C  THAT  -  SIGH  CB  TH22  HAS  XBSEBZED  DEC.  23,  1981 
C 

C  CALCULATE  AP20,  THE  POS1XI1E  TBAVELIIKG  BODE  AX 
C  ZoO. 

C 

UD  »  2.4Z|2)/(Z(1)*Z(2)) 

00  •  OU*(-OAT*K(l)*BETATH(l)*C11(1)*TB11) 

00  >  00/ |1AU(1)4JAT«K(1)) 

C 

TT  »  Tn22/(BETA*OAr*K(2)) tTH20/ (TA0(2) ♦0AI4K(2)) 
r»  ■  TV*JAI»K(2)  •cn(2)*eElAIH(2) 

TV  -  ((Z(1)-Z(2))/(Z{1)«Z(2)))*V¥ 

C  4-  . 

AP20  -  00  -  TV 
C 

C  CALCULATE  AIHO,  THE  HECATITE  TBATEILIRG  BOOB  AT 
C  THE  B0UNDAB1  Z-0. 

C 


I 


129. 

130. 

131. 

132. 

133. 

134. 

135. 

136. 

137. 

138. 

mo. 

141. 

142. 

143. 

144. 

145. 

146. 

147. 

148. 

149. 

150. 

151. 

152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 
161. 
162. 
163. 
163.1 

165. 

166. 

167. 

168. 

169. 

170. 

170.5 

171. 

171.5 

172. 
176. 
211. 
212. 

213. 

213.5 

214. 

215. 

216. 
220. 
221. 
222. 
226. 

229. 

230. 

231. 


00  •  TH23/(BE1A*JAY*K(2)) «1H20/(Tk0(2)  »JIt*K(2)) 

00  »  (2) ‘CIl  (2)  *EETATH  (2) 

00  -  00*  (-2.*Z(1)/(2(  1)  *Z(2))  ) 

■  C 

TT  «  (Z(1)-2(2))/(Z(l)*Z(2)) 

IT  ■  TV*JA1*1<  (1)  *BETATII(1)  *01 1  (1)  *181 1 
TT  »  VT/ (TAD (1)«JAY*K(1)) 

C 

AHIO  «  CO  *  TT 
C 

c 

C  MOV  CALCULATE  EAR  EIEID  BCDE  AHELITDCES  AP2E 
C  .  AMO  AflIP. 

C 

AP2P  »  JAY*K(2)*C11(2)*EE1A1H(2)*TB22 
AP2F  =  AP2E/(0AT*K(2)-BE1:A) 

AP2P  =  AP2F  ♦  JAY*K(2) *011(2) •BETATB{2)  *1820 
I  /(OAY*E(2)-TAO{2)) 

AP2F  a  AP2E  *  AP20 
C 
C 

AHIP  =  JAY*K  (1)*C11{1)  •BETATH(1)*IH11 
AnIP  =  AH1F/(JAY*K(1)-TA0(1)) 

AHIP  »  AHIP  *  AB10 
C 

C  CALCULATE  ACOUSTIC  POWER  AMD  STBAIS  IM  IBB 
C  PAB  FIELDS  CE  EEGIOKS  1  AMC  2. 

P2(ll)  =  CAES(AF2F*C0SJG  (AP2F) /8./ZC2J  ) 

P1(K)  =  CAES(AElE*C0HJG)AniE)/8./Z{1)) 

C 

C  .  • 

51  (MJ  *SOET  (2.  *2  (1) ‘PICS)) /Cl  1  (1) 

52  0I)  »  SCBT<2.*Z(2)*P2(A))/C11(2) 

C 

73  COMTINOE 
C 

C  .PRINT  OUT  SCEC  ANSVEBS 
C 

BRITE{6,367) 

367POBBAT('1  P  SI  .52 

1P1  P2  •/) 

DO  74  I  *  1,201 

NRITE(6.36e)  E(I),S1(I),S2(1)  ,P1(I),P2(1) 

74  CONTINUE 

368  FOBHATC  •  ,E  1P.B ,3I,E1  8. 8,3 X, El 8. B ,3Z, El8. 8, 3X, El 8.8) 

C 

C  CALL  PLOTTING  ROUTINES 
C 

CALL  STABTG( 'GENII, LI NEEX*8*», 0.0) 

CALL  SETSnG(23,1.0) 

CALL  SETSEG(24,1.0) 

CALL  SUPJEG(1.CEe,1.0E-20,1.0E10,1.0B>16) 

CALL  SETSF.C(132,1.) 

CALL  OBArHCCVEBIT*',  1,201, r,P2,0,'EBEQ**,»  I2*», 

1'1E-5  AIR///SI11CCN  AIEHA  =  3.29E6  12 

C  CALL  LJNISG{'VBRIT*',201,P,P)), 

CALL  EXITG 

C  VltAT  POLLCVS  IS  A  TABLE  CP  DATA  FOB  DIPTEBENT  BATEBIALS 
C  K  RFC  C  Ctl  BETATB 

C  PT8EX  DATA:  ADJUSTED  BETATB  1/8/80 


# 
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» 


232. 

C  .1260001  .2290Ct4 

.9750003  .73000011 

.270000-6 

233. 

C  HOLYBDENOn  tklK:  HOJUSTEO  BETATH  1/6/80 

234. 

C  .1790003  .1C22COS 

.1380003  .47000012 

.286000-5 

23S. 

C  COLD  DATA:  ACJDSTIC  EETATH  1/8/60 

236. 

C  .3450003  .1920CD5 

.1300003  .18900012 

. 130000-4 

237. 

C  WATER  DATA 

238. 

C  .6500000  .1000004 

;418S0D4  .23700010 

.690000-4 

239. 

C  ACETONE  DATA: 

240. 

C  .1980000  .7900CD3 

.2176004  .10600010 

.500000-3 

241. 

C  HETHANOL  DATA: 

242. 

C  .1550000  .7910003 

.2547004  .8300009 

.400000-3 

243. 

C  BELIUfi  (4K)  DATA: 

244. 

C  .200000-1  .  146CC03 

.52C0CD4  .7400007 

.1370000 

24S. 

C  ABGON  (85K|  CATA: 

246. 

C  . 1210000  .1417004 

.9840004  .10360010 

. 150000-2 

247. 

C  NITROGEN  (7SK)  CATA: 

248. 

C  .1380000  .8240003 

.1948004  .6380008 

. 170000-3 

249. 

C  SILICON  DATA:  ACJCSTED  BETATR  1/8/80 

250. 

C  .8350002  .23400C4 

.72800C3  .1657012 

.44800-5 

251. 

C  ALOHINUM  DATA:  AIOUSTED  BETATH  1/8/80 

252. 

C  .3000003  .2698904 

.9170003  .10900012 

.17800-4 

253. 

C  QUARTZ  DATA; 

254. 

C.  1350001  .220CC4 

.9370003  .86700011 

.55000-6 

255. 

C  GALLIUN  ARSENICE  CATA: 

256. 

C  .4600002  .530704 

.3500003  .11880012 

.59000-5 

257. 

C  CHR08IUS  DATA: 

258. 

C  .8700002  .7100C4 

.460003  .35000012 

.650000-5 

259. 

C  TANTALUH  DATA; 

260. 

C  .6800002  .  1660CS 

.142003  .26700012 

.650000-5 

261. 

C  SAPPHIRE  DATA:  ACOUSTED  BETATH  1/8/80 

262. 

C  .3300002  .39860C4 

.4330003  .49400012 

.300000-5 

263. 

C  CARBON  DATA:  AOJtSTED 

BETATH  1/8/80 

264. 

C  .8000002  .20000C4 

.7110C03  .32000011 

.520000-5 

264.1 

C  AIR  (GASEOOS)  CATA:  REV] 

264.2 

C  .2410000  . 1225CE1 

.7410003  .1416006 

.341000-2 

264.3 

C  IE-5  ATH  AIR  CATA  (CAICOLATED  PEOfl  1  ATB) : 

264.4 

C  .241000-5  .1225CD-4 

.7410003  .1416001 

.341000-2 

265. 

C 

a 

266. 

C  HHAT  EOLLONS  IS  A  TABLE  Or  ATTENOATICN  CONSTANTS  AND  SEPI.- 

267. 

C  ECTANCES  FOE  CERTAIN  RATEBJALS  1AIP  HCK,  SEC.  6G) 

268. 

C  BATERIAL  OPTICAL 

HAVE-  EOHEB  ATTEN  BEFLECTAKCE  - 

269. 

C  lENGTH, 

BICBCNS  CONST,  l/H 

(PCVEB) 

270. 

C 

271. 

C  ALfBVAP  .492 

1.4008 

.922 

272. 

C  AL,EVAP  .95 

1.  1208 

.912 

273. 

C  AO,ETAP  .50 

4.6007 

.504 

274. 

C  AO,eVAP  1.00 

7.6007 

.981 

275. 

C  SI,SHGL  CRST  .515 

2.4406 

.375 

276. 

C  SI, BULK  LC  EON  1.25 

2.9606 

.330 

277. 

C  S1,CVAP  .500 

1.8607 

.399 

278. 

C  5I,CVAP  .800 

3.2906 

.367 

279. 

C  HOLY,  BULK  .501 

7.5407 

.520 

280. 

C  BOLT,  BULK  1.000 

4.4607 

.574 

281. 

C 

321. 

STOP 

322. 

END 

* 

323. 

SDATA 

324. 

.10000E7  .10000E11 

1 

325. 

3.29E6 

326. 

.241000-5  .122500-4 

.7410003  .1416001 

.341000-2 

92 


93 


t 


a»»E«oxx  0 

PBOCBEB  TO  CALCOLATE  TBB  PHOTOACOOSTIC 
ttSPOUSZ  lU  A  3>SE(;iC»  SVSIEB. 


//DOPT3  JOO  •pn]E7«181*,*B  BBAT •  •CtASS«EtBEGIOH*S12X 
/•JOBPABfl  OESTsSEir 
//  EXEC  WATPIT 
//GO.SYSIM  DO  • 

SBATPIV  NOLIST 

C  THIS  POBTBIH  PBCGBAB  SOttlS  POB  TBE  ACODSTIC  STBAIB 
C  BAOIATEO  INTC  THC  SERI-INf INITE  HALF-SPACES,  BEGIONS 
C  1  AND  3,  AND  CENERAIED  BT  CPIICAL  AESCBBTICN  IN  BEOXOB 
C  2  WHICH  SEPABA1ES  TEEN  AND  HAS  THICKNESS  L.  THE 
C  INCOMING  OPTICAL  POKES  IS  ASSUMED  TO  BE  INCIDENT 
C  THROUGH  TBANSFABENT  HEOIOH  1.  AM C  IS  ASSOHED  TO  BATE 
C  OBIT  INTENSITY. 

C  HEDIUN  3  IS  ALSO  TBANSPABENT.  REIIDfl  t  EXTENDS  PBOB 
C  X«-IHPIKITT  TO  X=0.  HEDIDE  2  EXTENDS  PBOB  X=0.  TO  X>U 
C  BEDIDH  3  EXTENDS  FROM  X  =  L  TO  INFINITY.  THE  HEAT  PLOW 
C  IS  CALCULATED  USING  THE  BETHCD  CP  COUPLED  BODES 
C  AS  HELL.  THE  ACOUSTIC  GENEBATICH  IS  SOLVED  BY  TBE 
C  BETHOD  OP  CCDPIED  MOOES  IN  A  DISTfilBUTED  FINITE 
C  SOURCE  REGION  AS  PBESEKTED  IN  E. A.  AULD'S  TEXT  V.I,  CH.6. 
C  THIS  VEBSICN  CALCULATES  PCS  A  RANGE  CF  TFICRHESSES 
C  OP  REGION  2  FROM  0  TO  6  THERMAL  SKIN  DEPTHS. 

C  IT  THEN  PRODUCES  A  VERSATEK  PICT  CP  STRAIN  INTO  3 
C  VERSUS  A*L,  A  «  THERMAL  DIFFUSION  COEEFICIENT, 

C  L  s  LENGTH  OF  BEGION  2. 

C  THIS  VERSION  PICTS  POWER  INTC  BEGICB  3  P3  IBD 
C  POWER  INTO  BFGICN  1  PI.  A  TAFLF  OF 
C  DATA  IN  FORMAT  READY  TO  USE  (BEHOVE  COMHENT  C  IB 
C  COLUBN  1)  IS  STORED  IHHEDIATELY  BEFORE  THE  DATA  CARD. 

C  TO  BUN  OTHER  DATA,  CHANCE  THE  PLOT  LABELLING 
C  IN  THE  CALL  TC  GEAPKG;  CHANGE  THE  X  AND  Y  RANGE  IB 
C  TBE  CALL  TO  SUEJEG;  TRY  A  CHANGE  IN  THE  'LINEEX^S' 

C  OPTION  IN  THE  CALL  TO  STABTG  IP  AN  EBBOB  IN  TBB 
C  STARTG  CALL  RESULTS.  (E.G.  UIKEBlsfi') 

C 

REAL  SI  (2,200),  S3(2,20C),  STB(200)  ,  TADI(200), 

1  P1(2,200),P3(2,2C0),PHD(200) 

REAL*8  RHO(3),  C(3),  C11(3),  B£TATB(3) ,  ALPBi(3), 
1AA(3}  ,  20(3)  , 

2K(3),  KK(3),  L,LHAX,LINCB,CEL21,OEL23 
COBPLEX*16  TAU(3),  ABl 2, AP3I, AP32, 

1  JAY.  C1,C3,C20,C22, 

2  BETA,  2TH(3),  TH02, A , B ,CO,D2I ,D23,T81 1, 

3TH30,TU20,Tn2  1,TH22,U0, VV,HB,CDA 

C 

C  PRBQ  «  PBFCUENCY  OF  HODULATIOH  OP  LIGHT 

C  BET  -  OPTICAL  AHSORETICN  CCEPFXCIENT  OP  2 


READ,  FREQ 
HBITE(6,1IIO)PREO 
140  rOHBATCI  PBECOPHCT 
READ, BET 


«  ••K20.a,*  HZ') 


94 


56. 

57. 

50 

58. 

59. 

60. 

30 

61. 

62. 

63. 

C 

64. 

c 

65. 

c 

66. 

c 

67. 

c 

6B. 

c 

68.1 

c 

68.2 

c 

69. 

c 

70. 

71. 

72. 

40 

73. 

45 

74. 

75. 

76. 

77. 

47 

79. 

C  1 

80. 

c 

81. 

c 

■KITE (6, 50)  BET 

roKHiTco  onicAi  ios  coErr  bet 
BETA  -  CnilZItlT.O.) 

■R ITE(6,3C| 

roan AT ('0  kk 

1  C11 

00  «S  I«1,2 


E20.B,* 


BBTATB*) 


XK(I)  B  THEBEAt  CONDOCIXTITI  CF  ITH  BEGXOI 
CEMSJTY 

C(I)  B  SPECIFIC  BEAT 
C11(I)B  ELASTIC  STIFFNESS  C11 

BETATH(I)  B  lhERHAL  CCBFFICIEIT  CE  LINEAR  EZPAVSXOB 
■Oil:  AEJUSl  UANEBOOK  VALUE  CF  BEIATH 
BI  EUITIPLYING  El  EUIK  HCDOXUS,  DIVIDING  BY  C11. 

READ,  KK(I),  RHO(I).  C  (1)  ,  Cl  1  (I)  ,  BETATHd) 

■  RITE(6,<IC)  KK(I),  BUOil),  C(I),  CII(I},  DETATB  |I) 

FORSAT (•0',SE30.8) 

CONTINUE 
DO  180  N  B  1,2 

READ,  KK(i)  ,RRC(3)  ,C(3)  ,C11  (3)  ,EETATB(3) 

«R1TE(6,47)  KN(3),BHO(3),C(3),C11 (3),BETAIB(3) 

FORHATC  1',SE20.8) 

CALCULATE  AtPHA (I)  b  thEBEAL  EIFFCSIVITY  CF  XTR  LAYER 
AA(I)  B  tHEBBAL  diffusion  COEFFICIENT 
TAO(I)  =  C1»J)*AA(I) 

EO(T)  «  ACOUSTIC  IBPIEANCE 

K<i)  B  acoustic  propacmion  constant 

ZIR(I)  B  THEBBAL  ICPEDANCB 
02!  *  ZTa(2)/21B{1) 

023  B  ZIR(2)/Z1H(3) 

01121  «  Z0(2)/Z0i1) 

DEL23  B  Z0(2)/Z0C3) 

JAYb(0.,1.) 

DO  60  1=1,3 

ALPHA(I)  B  RR(I)/(RHOtI)*C(I)). 

AA(I)  =  DSQRT<(2.*3.  1<I15S*FPEC)/(2.*ALPBA(I))) 

ZO(I)  =  DSCRl  (BHC(I)*C11(I)) 

TAU(I)  B  CCnFLZ(AA(I1  ,AA(1)) 

K  (I)  =  2.*3- 14 159*FBEC*DSCBT (RHC (I) /C11 (1)  ) 

ZTH(I)  =  (1.*JAI)*0SQKT(3.14159*FBEC*BBC(I)*C(I)*KK(I)) 
■RITE(6,560)  1 
0  FORBATCO  I  B  ',  II) 

HRXTE(6,5e2)  AIPNA (I) ,AA (I) ,Z0  (I) ,T{I) 

2  FORBATC  ALFHA(I),AA(I),Z0(1),K(1)  =  i,  £018.6) 

HRITE(C,5E4)  T AU  (I) ,ZTH (I) 

4  FORBATC  TAO<I)  ,ZTH(I)  =  •,  2018. 6, 3Z, 2018. 6) 

VRITE(6,566)  ZTN(I) 

6  FORNATC  ZTH(I)  =  «,  2D18.6) 

CONTINUE 


ZTH(I) 


021  B  ZTH|2)/ZTB(1) 
023  B  ZTH(2)/ZTH(3) 
0EL21  B  Z0(3)/Z0<1) 
0EL23  B  ZC(2)/Z0<3) 
HRITE(6,£68)  121,023 
FORBATC  021,023  =  t 


«BITE(6,570)  CEL21,0EL2J 


2016.6,31,2018.6) 
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IH. 

117. 

118. 

119. 

120. 
122. 
123. 
12«. 

125. 

126. 

127. 

128. 

129. 

130. 

131. 

133. 

134. 

135. 

136. 

137. 

138. 

139. 

140. 

141. 

142. 

143. 

144. 

145. 

146. 

147. 

148. 

149. 

150. 

151. 

152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 
161. 
162. 

163. 

164. 

165. 

166. 
167. 

169. 

170. 

171. 

172. 

173. 

174. 

175. 

176. 

177. 


570  rOIBATf*  CBL21,ttL23  «  1018.6,31,1018.6) 

c 

TSO  •  1./44I2) 

HBITe|6,6<)  ISO 

65  rOBRATCO  TfllBflAL  SKIN  tttXH  1/AA  (2)  «  *,  820.8} 

C 

C  SET  OP  LOOP  TO  IHCOEREHT  LtHGTE  7008  0  TO  IHil, 

C  HHERB  LBAK  IS  SUCH  THAT  k(2)*IBAX  »  6.,  IH  200  STEPS 
C  BECALL  1/AA  IS  THE  THEBHAL  SKIH  DEPTH. 

C 

CALL  TBAPSfl, 1,10000,1,1) 

LHAX  «  6.0/kA(2) 

LZVCB  «  LHAZ/300. 

L  -  0. 

DO  200  8  >  1,  200 

L  •  L  ♦  IIHCB  - 

C 

C  CALCULATE  TRRPERATUBE  AEPIITODES  TH11,  TH30  IB  BEGI08S 
C  1  AHO  3,  AND  THEK  TH20  (EXP (~T AD«Z) ,  TB2 1  (EXP (tTAD*Z) ,  ABD 
C  TH22 (EXP (-BETA*Z)).  TH02  IS  A  UBICDITOUS  COLLECTION  OF 
C  CONSTANTS,  AS  ARE  A,  B,  AKC  CD  (==  COHECN  CEKOHINATOB). 

C 

TB02  *  B£7A/(TA0(2)*ZTH(2)) 

T802  -  TBC2/((EETA*BETA)/(TAD(2)*TAD(2))*1.) 

C 

A  «  BBTA/(BETA«TkD(2)) 

A  »  A*(CCEXP(-(BETA*TA0(2))  *L)-1,) 

C 

B  •  BBtA/(£ETA-TA0{2)) 

B  «  B*(CDEXP(-BETA*L)  -  CDEXP  |-TAtJ  (2)  •!)  ) 

C 

CD  »  (1.-D21) • (1.-D23)*CDEXP{-TI0(2)*1) 

CD  *  CO  *  {1.«C21)«(1.*C23)*CDEXP(TA0{2)«t) 

C 

TH11  »  B*{1.-D23)  -  A*(1.*D23)»CDEXI(TA0C2)*L) 

TH11  =  -TB 11/ (218(1)  ♦CO) 

C  -  . 

TH30  =  D*(1.*C21)*CDEXE(TAD(2)*I)  -  A*  (1.-021) 

TH30  >  TH30/(2TB(3)*CO) 

C 

C 

TH20  -  (1.-D23)*(1.-(D23*BETA)/TA0(2))*CDEXP(-BETA*L) 

TH20  >  TH20  -  (1. «D23) •  (1. «  (D21*BE1 A) /TAO (2) ) •CDEXP (TAD (2) *1) 
TH20  >  TR20«(TE02/C0) 

C 

TH21  »  (1.-C23)*(1.*(D21*BETA)/TA0(2))*CDEXP(-TA0(2)^L) 
T821«TH21-  (1.*D21)*(1.-(D23*BE1A)  /IAU(2)  )  •CDEXP  (-6£TA*L) 

TH21  ♦  TB21*  (TH02/CD) 

C 

TH22  -  -TH02 

C  -  - 

c 

C  CALCULATE  Cl,  C20,  C22,  AND  C3,  WHICH  ABE  COHHOH 
C  COLLECTIONS  CP  CCNSTANTS  SEEN  IN  THESE  PABTS. 

C 

Cl  ■  -JAI*J'(1)*C11(1)*BETATH(1)*TB11 
Cl  ■  Cl/(JAr*IC(1)  ♦  TAO(I)) 

C 

C3  •  -JkY*K(3)*Cl1(3)*eETATB(3)*TH30 
C3  •  C3/(JAT*A(3)  ♦  TAO  (3)) 
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178. 

179. 

180. 
181. 
182. 
183. 
18«. 
IBS. 
186. 

187. 

188. 

189. 

190. 

191. 

192. 

193. 
199. 

195. 

196. 

197. 

198. 

199. 

200. 
201. 
202. 
203. 

209. 

205. 

206. 

207. 

208. 

210. 
211. 
212. 
213. 
219. 

215. 

216. 

217. 

218. 

219. 

220. 
221. 
222. 
223. 
229. 

225. 

226. 

227. 

228. 

229. 

230. 

231. 

232. 

233. 
239. 

235. 

236. 

237. 

238. 


C 

00  •  CDeXP|-(J>Y*8(2)»T9af2))*l)  -  1. 

00  •  aO*(1H20/(JAY«K(2)*190(2))) 

'C 

TT  «  CDEIP(-(0AT*f  (2J-T«0<2))  *11  -  1. 

ft  -  V«*(1M31/(JAt*X(2)-1iOC2n) 

C 

00  «  CDEXF{-(CAt«K(2)»BETA)*t}  -  I. 

00  «  0B«(1H22/(JAt*X(2)*eEIA)) 

C 

C20*3AT*K  (21 *011 (2|  *B£TA1H (2) * (D0»Vf*S8) 

C 

c . 

00  ■  CDEXP((3AX*l<{2)-TA0(2))*l)  -  1. 

00  «  00«(TU20/(JAX«E(2)-TA0|2})) 

C 

fT  »  CDEXP((3Al*Ji(2)  ♦TAD(2))*t)  -  1. 

TX  e  VV*(TE2V(JAX*K(2HTAa(2))) 

C 

00  *  CDEXP  ( (JAY*P(2)-BEtA)  *11  -  I.  i 

00  «  0M*(Iil22/|JAY*X(2)-BElA)) 

C 

C22  =  -3AY*B (3)*C11 (2)  *EE1A1H(2)*CDEXP(-JAY*K(2)*L) 

C22  «  C22*(00*VV«0B) 

C 

C  THERE  ARE  TOC  COEttOH  DEACfllSAICBS.  THY  THBBBAt 
C  ORE  HE  CAILIt  CC.  THE  ACOOSTXC  OHE  HE  SHALL  CALL  CDA 
C 

COA  «  (1.4DE121J*(1.+DIL231*C0EXP (JAT*K(2|*l) 

COA  «  (1.-t)E121)*(1.-DEl23)  •CBEXP(-JAI*K(2)91)  -  CDA  . 

C 

C 

C  CALCULATE  AH12,  THE  PAB-rXElO  BCDE  MBPlXTODg  II 
C  BEDXOB  1 
C 

00  *  -C1.-CE123)*C22  -  2.*CIL23*C3 

00  ■  -  (1. ♦DE123J *C20*CDEXP(JAY*P(2) ♦L)  *00 

TV  »  (1.4CEI23)*(JAr»K(1)  ♦  DEL21  •!  AO  (1)  )  *CDEXP  (JAY*K  (2)  *1) 

00  »  (1.-CEL23)*{JAY*KC1X  -  D£L21*7AU  (  1)  )  *CDEZP  (-JAX*K  (2)  *6) 

00  *  00  -  TV 

00  ■  -C1*00 

00  »  00/(OAT*F(1)-TA0{1)) 

C 

AHU  «  00  4  00 

AH1Z  «  AHIZ*  (2./C0A)  , 

C 

c 

c  BOH  SOLVE  POR  AP3t.  THE  YALOE  CP  THE  PCSITITE  TBAVELLXOG 
C  BODE  XB  HEClUB  3  AT  THE  BCUNDABY  OITH  BEDIOH  2 
C 

00  >  .5*(1.4DEL21)*(2.*C22-(1.-tEL23)*C3)4CDBZP(JAZ*K(2l*'L) 
VV  ■  .5*  (1.-DE121)4(1.40EL23)  •C34CCYXP(-JAY*KC2J*L) 

DO  ■  2. 407121*01  4  (1.-0ei2l}4C20 
C 

AP3t  ■  -2.*((CO*VV*OHI/CBAi 
C 

C  OOO  CALCULATE  APJZ,  THE  PAB  PXEtl  RCCE  ABPLXTDDE  Xl 
C  HEDXOH  3  . 

C  4 

AP3Z  •  JAI*E(3i4Cl1(3}  •BETA1H(3}41B30 


239.  AV3Z  -  A»3Z/(JftT*K(3)-t«0f3))  •  ZP3L 

290.  C 

263.  C 

269.  C 

266.  C  CtLCOLITl  STBIIM  St  SkCliTfC  1910  RI6Z0I  1,  AID 

267.  C  SI  lADIATEO  lilTC  REGION  3.  DRIT  ItTZIISITl 

268.  C  IS  ASSOHZC  fOB  TEE  IBCIOEBI  LICBl. 

269.  C 

270.  202  >  SB6t(CDABS(AniZ)/(2*cn(1))) 

271.  S3(N,H)  «  SBCI(CDABS(AE3Z)/(2*C11(3))) 

272.  PI(N.H)  >  SNGl(CDADS(ARtZ*DCOKJG(AP.lZ)/8./Z0(1))) 

273.  P3(M,n)  c  SHG1(CCA£S(AP3Z*CC0NJC(AP3Z)/E./Z0(3||) 

279.  C 

275.  C  LOOP  ONTIL  200  IRCBZHEBIS  Of  LBRCIB  BABE  BBEB  BADE 

276.  C 

277.  TAOLCH)  «  SBGt (AA (2)  •!} 

278.  200  COKtIHUE 

279.  NRItE(6,1S0| 

280.  150  POBHATCO  PCEEB  IBTO  1  ECBBB  INTO  3*> 

281.  00  170  fl«1,20C 

282.  HRITE(6,  UO)  F1(K.B)«  P3(B,B) 

283.  160  rOBBATC  •«  2120.8} 

289.  170  COKTIHUE 

285.  180  COBTIHOE 

286.  C  . 

287.  C  CALL  PLOTTIBG  BOOTIRES 

288.  C 

289.  CALL  STAB1G(*CEN11.LIHBBZ«6**.0.0) 

290.  CALL  SETSeG(29.1.0) 

291.  CALL  $U8JEG(0.0. 1.0E-20, 6.00,1. OE'IS) 

292.  CALL  SETSEG  (132,1.) 

293.  DO  220  3«1,20C 

299.  PHR(3)  *  F3(1,J) 

295.  220  COKTIHUE 

296.  CALL  GRAPEG(*TERIT**,1,200,1ADL,PHB,0,'AI**,'X2,  B/B**, 

297.  I'SAPPHIBE/SIIICCH  //  HA1EB  13  1  GHZ*') 

298.  00  230  J  >1,2C0 

299.  PVR(J)  »  E3(2,3) 

300.  230  COKTIHUE 

301.  CALL  LINESC('TERIT«S200,TADL,PH) 

302.  CALL  EXI1G 

30«.  C  HRAT  rOLLCVS  IS  A  TABLE  Of  EATA  FOR  OlFEEBEHl  BATEBZALS 


305. 

C 

K 

BBC 

C 

C11 

EETATB 

306. 

c 

PZBEX  DATA: 

ACJUSTED 

BETATB  1/8/80 

307. 

c 

.1260001 

.2290CD4 

.9750003 

.73000011 

.270000-6 

308. 

c 

BOLTUDENUB 

DAK:  AtJOSTED  EETAIE 

1/8/80 

309. 

c 

.1790003 

. 1022005 

.1380003 

.47000012 

.286000-5 

310. 

c 

COLO  DATA: 

AtJCSTIt  DETATH  1/8/80 

311. 

c 

.3450003 

.1920005 

.1300003 

.18900012 

.130000-9 

312. 

c 

BATER  DATA 

313. 

c 

.6500000 

. 1000004 

.9185009 

.23700010 

.690000-9 

319. 

c 

ACETONE  DATA: 

315. 

c 

.1980000 

.7900003 

.2176004 

.10600010 

.500000-3 

316. 

c 

BETHANOL  DATA: 

. 

317. 

c 

.1550000 

.7910CD3 

.2597004 

.8300009 

.900000-3 

318. 

c 

BELIOB  (4K) 

CATA: 

6 

319. 

c 

.200000*1 

.1460003 

.5200009 

.  .7400007 

.1370000 

320. 

c 

ARGON  (85K) 

DATA: 

321. 

c 

.1210000 

. 1417CD4 

.9890009 

.10360010 

.150000-2 

322.  C  HZTROGEN  (75R)  EATA: 


150000*2 


323. 

C 

.1380000 

8240C03 

.1948004  .6380008 

.170000-3 

32«. 

c 

SILXCOH  04T4: 

ADJUSTED  BETATH  1/8/80 

325. 

c 

.8350002 

.2340004 

.72H00C3  .1657012 

.44800-5 

32«. 

c 

ALOniHUN  0411 

i:  AtJOSTED  BETATH  1/B/80 

327. 

c 

.3000003 

.2698904 

.91700C3  .10900012  .17800-4 

328. 

c 

QOtBTZ  DATA: 

329. 

c 

.1350001 

.220004 

.9370003  .86700011  .55000-6 

330. 

c 

CALLIOH  ARSCMICE  tATA: 

331. 

c 

.4600002 

.S307C4 

.3500003  .11880012  .59000-5 

332. 

c 

CHBOHIOn  OAT A 

• 

• 

333. 

c 

.8700002 

710014 

.460CD3  .35000012 

.650000-5 

334. 

c 

TAVTALOfl  DATA 

335. 

c 

.6800002 

166QE5 

.142003  .26700012 

.650000-5 

336. 

c 

SAPPHIRE  DATA 

:  AEJUSTEO  BETATH  1/8/80 

337. 

c 

.3300002 

39860C4 

.4330003  .49400012 

.300000-5 

338. 

c 

CABBOH  DATA: 

ACJCSTED 

1  BETATH  1/8/80 

339. 

c 

.8000002 

2000CC4 

.7110003  .32000011 

.520000-5 

340. 

c 

AIB  (GASEOOS) 

DATA:  ME8I 

341. 

c 

.2410000 

12250E1 

.7410003  .1416006 

.341000-2 

342. 

c 

IE-5  ATB  AIR 

EATA  (CALCULATED  PBOB  1  ATB) : 

343. 

c 

.241000-5 

122500-4 

.7410003  .1416001 

.341000-2 

344. 

c 

345. 

c 

VHAT  P0LL08S 

IS  A  TABLE  OP  ATTEKDATIC8  COHSTAHTS  ABO  BEPL- 

346. 

c 

ECTAHCES  POP 

CIRTflH  HATE8IILS  (IIP  HBK,  SEC. 

6G) 

347. 

c 

BATERIAL 

CPT1CAI 

SAVE-  POBEB  ATTEB 

BEPLECTAHCB 

348. 

c 

ZEBGTU, 

BICBCSS  CCttSl,  1/B 

(PC8EB) 

349. 

c 

350. 

c 

AI.,ETAP 

.492 

1.4008 

.922 

351. 

c 

AL.BVAP 

.95 

1.12D8 

.912  . 

352. 

c 

AOfEVAP 

.50 

4.6007 

.504 

353. 

c 

A0,EVAP 

1.00 

7.6007 

.981 

354. 

c 

SI,SNGL  CRS1 

.515 

2.4406 

.375 

355. 

c 

SI,BULK  LO  FOB  1.25 

2.9606 

.330 

356. 

c 

SIfETAP 

.500 

1.8607 

.399 

357. 

c 

SI,ETAP 

.800 

3.2906 

•  367 

358. 

c 

BOLT,  BULB 

.SCI 

7.5407 

.520 

359. 

c 

BOLT,  BULK 

1.000 

4.4607 

.574 

360. 

c 

400. 

STOP 

t 

401. 

END 

• 

402. 

SDATA 

403. 

.10000010 

404. 

3.29E6 

405. 

.3300002  .: 

3886CD4 

.4330003  .49400012 

.300000-5 

406. 

.8350002 

.2340004 

.7280003  .1657012 

.44800-5 

407. 

• 

6500000  .1000004 

.4185004  .23700010 

.690000-4 

408. 

• 

6500000  .10000C4 

.4185004  .23700010 

.690000-4 

409. 

SSTOP 

99 
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